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ABSTRACT 


This  report  covers  in  detail  the  solid  state  research  work 
of  the  Solid  State  Division  at  Lincoln  Laboratory  for  the 
period  1  February  through  30  April  1970.  The  topics 
covered  are  Solid  State  Device  Research,  Materials  Re¬ 
search,  Physics  of  Solids,  and  Microelectronics.  The 
Microsound  work  is  sponsored  by  ABMDA  and  is  reported 
under  that  program. 
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INTRODUCTION 


I.  SOLID  STATE  DEVICE  RESEARCH 

Three  high  purity  epitaxial  GaAs  samples  have  been  evaluated  as  far  infrared  detectors  under 
different  background  conditions.  Spectral  response  measurements  at  4.  2°  and  1.93°K  indicate 
that  the  long  wavelength  response  from  200  to  350  p  is  due  primarily  to  a  combined  photo- 
thermal  excitation  process.  The  conductivity  and  thus  the  performance  characteristics  of  the 
detectors  consequently  depend  critically  on  both  the  detector  temperature  and  background 
radiation.  The  best  long  wavelength  response  was  observed  at  4. 2°K  in  reduced  background. 

Long  wavelength  laser  emission  at  28  p  has  been  observed  from  Schottky  barrier  diodes  fab¬ 
ricated  by  evaporating  Pb  contacts  on  Bridgman-grown  Pb^  ^Sn^  285Te*  ^le  materiai  usecl 
was  4.  5  x  1017cm-3  p-type  and  was  obtained  as  part  of  a  study  of  the  preparation  and  charac¬ 
terization  of  metal  saturated  Pbj  xSnxTe  as  a  function  of  isothermal  annealing  temperatures. 

N-p  junction  photovoltaic  detectors  have  been  fabricated  in  p-type  Hg^  ^-Cd^  ^Te  using  proton 
bombardment  to  create  the  n-type  region.  The  photovoltaic  response  of  these  diodes  peaks 
at  1.8pwith  a  quantum  efficiency  of  15  percent.  The  300°K  detectivity  at  the  peak  is 
2  x  1010cm  HzV2/w. 

A  vidicon  type  television  camera  tube  has  been  constructed  employing  an  array  of  reverse 
biased  GaAs  Schottky  barrier  photodiodes  to  form  the  image  sensing  target.  This  target  has 
many  advantages  over  a  conventional  photoconductive  vidicon  target,  and  its  simplicity  of 
fabrication  and  stability  make  it  an  attractive  alternative  to  the  silicon  p-n  photodiode  array 
target. 

A  new  technique  has  been  developed  for  the  fabrication  of  semiconductor  devices  utilizing  a 
thin  semiconductor  layer.  The  advantage  of  this  is  that  different  parts  of  the  devices  can  be 
fabricated  on  both  sides  of  the  layer  in  close  proximity  to  each  other  adding  an  additional 
dimension  to  the  fabrication  technology.  The  technique  has  been  used  to  fabricate  GaAs 
Schottky  barrier  field  effect  transistors  in  which  the  source  and  drain  ohmic  contacts  sepa¬ 
rated  by  5  p  are  on  one  side  of  an  0. 5- p  thick  layer  of  GaAs  and  the  5  p  long  Schottky  barrier 
gate  is  on  the  opposite  side  positioned  over  the  source-drain  gap.  This  structure  and  others 
utilizing  this  technique  should  result  in  improved  high  frequency  performance  of  many  semi¬ 
conductor  devices. 

Acousto-electric  amplification  coefficients  have  been  measured  at  room  temperature  over  the 
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frequency  range  from  1.23  to  3.75  GHz  in  5  x  10  cm  n-type  GaAs  using  Brillouin  scat¬ 
tering.  Good  agreement  with  theory  is  observed  for  the  absolute  value  of  the  gain  coefficient 
over  this  range  and  for  the  frequency  of  maximum  gain. 
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II.  MATERIALS  RESEARCH 

A  new  seed  rod,  cooled  by  means  of  a  heat  pipe,  has  been  used  for  pulling  crystals  in  the  tri¬ 
arc  Czochralski  furnace.  Use  of  the  heat  pipe,  a  device  in  which  extremely  efficient  heat 
transfer  is  accomplished  by  the  evaporation  and  condensation  of  a  volatile  liquid,  eliminates 
the  disadvantages  encountered  with  a  water-cooled  seed  rod. 

Single  crystals  of  CdO  with  dimensions  up  to  about  5x5x8  mm  have  been  grown  from  the 
vapor  phase  at  about  1050°C  by  an  open  tube  method  in  which  a  stream  of  helium  saturated 
with  cadmium  vapor  is  mixed  with  an  oxygen  stream  inside  a  cylindrical  liner  made  of  high 
purity  A^Og. 

Charts  of  the  standard  free  energy  of  formation  (AG^°)  as  a  function  of  temperature  have  been 
constructed  for  binary  selenides  and  tellurides  on  the  basis  of  published  thermodynamic  data. 
No  charts  of  this  type  have  previously  been  available  for  these  compounds. 

The  effects  of  deviations  from  stoichiometry  on  the  lattice  parameter  and  electrical  resistivity 
of  EuO  have  been  confirmed  by  measurements  on  europium-saturated  and  oxygen- saturated 
single  crystals.  A  marked  increase  in  the  optical  transmission  of  EuS  and  EuSe  single  crys¬ 
tals  grown  from  europium-rich  solutions  has  been  achieved  by  annealing  the  crystals  in  sealed 
tungsten  crucibles. 

Optical  transmission  measurements  on  ZnTe^_xSex  single  crystals  have  confirmed  the  exist¬ 
ence  of  an  energy  gap  minimum  in  this  alloy  system.  Similar  measurements  demonstrate  a 
strongly  nonlinear  dependence  of  energy  gap  on  composition  in  the  ZnTe^  ^Sx  alloys,  although 
no  minimum  can  be  observed  because  the  mutual  solubility  of  ZnTe  and  ZnS  is  too  low. 

Optical  transmission  measurements  on  CdTe^_xSex  single  crystals  have  established  the  exist¬ 
ence  of  an  energy  gap  minimum  in  this  alloy  system.  The  data  indicate  that  the  energy  gap 
for  a  given  composition  is  about  0.03eV  higher  for  samples  with  the  hexagonal  wurtzite 
structure  than  for  those  with  the  cubic  zincblende  structure. 

An  energy  level  scheme  which  takes  account  of  the  role  of  cation  vacancies  has  been  proposed 
for  NiO.  This  scheme  accounts  for  the  optical  and  electrical  properties  of  both  undoped  and 
Li-doped  NiO. 

The  polaron  morphology  in  vanadium  spinels  has  been  investigated  by  analyzing  the  observed 

9+ 

variation  in  Seebeck  coefficient  with  composition  in  the  systems  Mg[Bx  xlO^  where 

B  =  Mg  or  Ni.  The  results  indicate  that  with  decreasing  V-V  separation  (R),  the  polarons 
44- 

change  from  V  cations  to  X^-cluster  and  V^,-cluster  polarons,  and  finally  to  itinerant  elec¬ 
trons,  in  a  narrow  range  2.99  >  R  >  2.94  A. 

Integrated  intensity  data  obtained  in  x-ray  diffraction  measurements  on  powdered  samples 
have  been  used  to  carry  out  a  structure  refinement  for  the  hexagonal  (6H)  form  of  RbNiF^. 
The  results  are  generally  in  agreement  with  published  values,  but  there  are  significant  dis¬ 
crepancies  in  the  positions  of  the  fluorine  atoms. 

The  effect  of  high  pressure  on  the  crystal  structure  of  CsMnCl^,  RbMnCl^,  and  RbFeCl^  has 
been  determined  by  means  of  x-ray  diffraction  measurements.  For  each  compound,  the 
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proportion  of  hexagonal  close  packing  decreases  with  increasing  pressure,  and  the  ultimate 
high  pressure  form  has  the  perovskite  structure,  which  is  entirely  cubic  close  packed. 

The  sources  of  error  in  measuring  volume  compression  by  means  of  x-ray  diffraction  studies 
at  high  pressures  have  been  analyzed.  The  procedures  required  to  obtain  accurate  lattice 
parameters  for  the  sample  and  for  a  marker  material  such  as  NaCl  used  to  determine  the 
applied  pressure,  are  discussed  in  detail. 

The  volume  compression  of  TiOx  for  samples  with  x  between  0.83  and  1.25  has  been  inves¬ 
tigated  by  means  of  x-ray  diffraction  studies  made  with  diamond-anvil  high  pressure  units. 

-4 

For  pressures  up  to  140  kbars,  the  measured  values  of  — V/VqAP  are  between  2.3  x  10  and 
4.0  x  10-4  kbar-1. 

III.  PHYSICS  OF  SOLIDS 

The  quasi-particle  approach  to  the  optical  properties  of  simple  metals  has  now  been  applied 
to  the  calculation  of  the  optical  properties  of  potassium.  Using  nonlocal  pseudopotentials  and 
the  experimental  Fermi  surface,  both  the  photon  energy  threshold  for  inter  band  absorption 
and  the  magnitude  of  the  absorption  have  been  calculated. 

Work  on  T^O^  continues.  In  an  attempt  to  describe  the  nature  of  the  insulator-metal  transi¬ 
tion,  a  model  based  on  the  competition  between  elastic  and  electronic  free  energies  of  the 
system  as  a  function  of  temperature,  is  being  investigated.  Although  this  model  is  capable  of 
giving  a  reasonable  description  of  the  change  in  c/a  ratio  with  temperature,  it  predicts  a  very 
pronounced  "mode  softening”  in  the  vicinity  of  the  transition  which  is  inconsistent  with  our 
Raman  scattering  data. 

A  mechanism  for  the  description  of  a  metal  to  semiconductor  transition  with  increasing  tem¬ 
perature  has  been  found.  The  result  is  based  on  a  comparison  of  the  free  energy  for  localized 
(Wannier)  vs  band  (Bloch)  states  for  a  half-filled  band  Hubbard  Hamiltonian. 

The  phonon  pinning  spectrum  of  an  electron  bound  to  an  impurity  in  a  semiconductor,  as 
affected  by  the  presence  of  the  impurity,  is  under  study.  A  preliminary  calculation  indicates 
pinning  due  to  interaction  with  local  or  resonant  modes,  as  well  as  pinning  due  to  lattice 
modes. 

A  band  theory  of  optical  activity,  which  arises  from  terms  in  the  dielectric  constant  linear  in 
q,  has  been  developed.  The  general  results,  when  reduced  to  a  simple  four-band  model,  give 
a  value  for  the  ratio  of  rotary  power  parallel  and  perpendicular  to  the  c-axis  which  is  close  to 
our  previously  measured  ratio  in  o-quartz. 

Work  on  the  magnetic  semiconductor  EuO  has  now  been  extended  to  studies  of  magnetic  prop¬ 
erties  in  the  vicinity  of  the  Curie  temperature.  Measurements  have  been  analyzed  to  deter¬ 
mine  the  magnetic  critical  indices. 

In  other  magnetic  work  the  study  of  two-magnon  Raman  scattering  in  KNiF„  has  been  extended 

2+  6 
to  samples  doped  with  Mn  impurities.  In  addition  to  the  intrinsic  two-magnon  excitation 

previously  observed  in  pure  samples,  anew,  relatively  narrow  line  has  been  observed,  with 
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an  intensity  that  is  dependent  on  the  Mn  concentration.  The  line  is  tentatively  identified  as 

2+ 

a  resonant  two-magnon  state  associated  with  an  Mn  impurity.  The  Raman  spectrum  of  the 

o 

antiferromagnet  CsMnF^  has  been  observed  using  the  4880-A  radiation  from  an  argon  laser. 
In  addition  to  the  active  phonon  modes,  a  mode  has  been  observed  whose  shape,  position, 
temperature  dependence  and  polarization  indicate  that  it  arises  from  the  excitation  of  two 
magnons. 

Theoretical  work  on  light  scattering  from  acoustic  plasma  waves  and  single  electron  excita¬ 
tions  in  magnetoplasmas  continues.  A  model  calculation,  involving  only  scattering  due  to  den¬ 
sity  fluctuations  and  ignoring  spin-orbit  coupling,  has  been  carried  out  to  illustrate  some  of 
the  general  features  of  magnetoplasma  quasi-elastic  scattering. 

In  another  phase  of  the  laser  program,  heterodyne  experiments  with  both  millimeter  and 
infrared  sources  have  been  performed  in  high  purity,  epitaxial  GaAs.  The  millimeter  waves 
produce  avalanche  breakdown,  and  mixing  via  hot  electrons  is  observed.  At  laser  frequencies 
the  DC  bias  level  determines  whether  detection  proceeds  via  a  similar  hot  electron  mechanism 
or  through photothermal  ionization  of  shallow  donor  states.  The  time  constant  associated  with 
the  latter  process  is  significantly  faster  than  that  of  hot  carrier  relaxation. 

Two  nonlinear  effects  on  the  polarization  and  linewidth  of  an  intense  broadband  maser  signal 
have  been  calculated  with  application  to  the  interstellar  masers.  One  is  the  competition 
between  opposite  circularly  polarized  Zeeman  modes  by  means  of  parametric  down-conversion 
involving  an  electron  cyclotron  idler  wave.  The  second  is  the  saturation  broadening  in  a  two- 
level  system  due  to  a  broadband  signal  having  Gaussian  statistics. 

IV.  MICROELECTRONICS 

The  Microelectronics  Group,  which  is  now  part  of  the  Solid  State  Division,  has  a  dual  charter 
of  service  and  research  in  integrated  electronics.  Service  tasks  have  provided  an  opportunity 
to  apply  new  processes  and  techniques  that  have  been  developed  as  part  of  the  research  pro¬ 
gram  on  integrated  circuits.  Typically,  at  any  one  time  there  are  about  20  active  programs 
ranging  from  special  semiconductor  devices  to  large  hybrid  integrated  assemblies.  In  addi¬ 
tion,  about  five  or  six  research  programs  are  in  progress  as  long  range  support  for  Lincoln 
Laboratory  programs. 

The  present  programs  in  research  and  development  include  the  following: 

(1)  Rigid  and  nonrigid  forms  of  beam  lead  substrates  have  been  fabricated  with 
gold  and  more  recently  with  aluminum  metallization  for  the  beams  and  air 
gap  crossovers.  Some  work  remains  to  be  done  to  improve  yield  and  cost 
reduction  in  the  fabrication  and  assembly  of  the  substrates,  but  the  basic 
techniques  are  developed  for  several  forms  of  substrates  including  an 
automated  version. 

(2)  Aluminum  beam  leaded  silicon  chips  have  been  developed  as  a  comple¬ 
mentary  assembly  method  to  the  beam  lead  substrates.  Air  gap  cross¬ 
overs  using  aluminum  metallization  on  the  silicon  chips  also  look  prom¬ 
ising  as  a  means  of  providing  an  additional  and  high  yield  interconnection 
method  at  the  chip  level. 
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(3)  A  silicon  dioxide  passivation  method  for  mesa  semiconductor  devices  has 
been  developed  which  has  several  advantages  over  present  passivation 
methods.  Actual  devices  have  been  fabricated  using  this  technique  with 
excellent  results. 

(4)  A  development  effort  on  Impatt  and  Trapatt  microwave  devices  with 
special  geometries  has  been  pursued  during  the  past  year.  New  packaging 
and  die  attachment  techniques  are  being  investigated  to  reduce  the  thermal 
coefficient  and  improve  the  efficiency. 

(5)  A  laser  scanner  for  mask  inspection  and  semiconductor  chip  evaluation 
has  been  built  and  tested.  Improvements  in  the  scanning  mechanism  are 
under  way  which  will  allow  faster  scanning.  Metallization  faults ,  mask 
misalignment,  diffusion  faults  and  oxide  problems  and  device  gain  char¬ 
acteristics  can  be  examined  with  the  scanner. 

Other  programs  that  are  more  applications  oriented  and  include  some  research  or  develop¬ 
ment  are  listed  below. 

(1)  A  strip  or  an  array  of  64  silicon  photodiodes  spaced  6  mils  on  centers  has 
been  developed  for  a  satellite  stabilization  system.  Two  approaches  have 
been  taken:  one  method  employs  4  strips  of  64  diodes  assembled  very 
accurately  in  a  single  package,  and  a  second  method  employs  a  redun¬ 
dancy  and  discretionary  wiring  to  allow  the  fabrication  of  the  four  diode 
arrays  on  a  single  silicon  wafer. 

(2)  Special  diodes  with  very  shallow  junctions  (0. 5  ji)  have  been  fabricated  for 
use  in  high  power  applications  with  an  electron  beam  as  the  input  control. 

These  devices  hopefully  will  compete  favorably  with  magnetrons  and 
traveling  wave  tubes  in  power  and  frequency.  Present  results  indicate  an 
oxide  charging  problem  and  lower  power  capabilities  than  anticipated. 

(3)  A  multi -element  particle  detector  using  a  silicon  array  is  under  develop¬ 
ment  in  cooperation  with  Dr.  Kosta  M.  Tsipis  (Laboratory  for  Nuclear 
Science,  M.I.T.). 

Other  programs  of  a  service  nature  that  require  moderate  development  are: 

(1)  Surface  wave  transducer  program. 

(a)  Fine  line  masks. 

(b)  Semiconductor  amplifier  substrates. 

(2)  Microwave  thin  film  assemblies. 

(3)  Mixer  diodes  employing  5-jiSchottky  barrier  devices. 
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I.  SOLID  STATE  DEVICE  RESEARCH 


A.  GaAs  FAR  INFRARED  DETECTORS 

Three  high  purity  epitaxial  GaAs  samples  have  been  evaluated  as  far  infrared  detectors 
under  two  different  background  conditions.  The  material  properties  of  these  samples  are  given 
in  Table  1-1.  Although  the  samples  fall  in  the  range  of  purity  generally  obtainable  today  by  epi¬ 
taxial  techniques,  two  of  them  have  relatively  low  compensation  ratios  compared  to  those  usu¬ 
ally  observed. 


MATERIAL 

TABLE  1-1 

PROPERTIES  OF  DETECTOR  SAMPLES 

nd 

na 

n300°K 

h77°K 

En 

Sample 

14  -3 

14  -3 

14  -3, 

,  2  . 

D 

No. 

(10  cm  ) 

(10  cm  ) 

(10  cm  ) 

(cm  /V-sec) 

(meV) 

1 

1.2 

1.0 

0.  13 

155,000 

5.4 

2 

2.0 

0.4 

1.6 

153,000 

5.  1 

3 

4.9 

1.0 

3.7 

108,000 

4.6 

The  experimental  arrangement  shown  in  Fig.  1-1  was  used  for  the  optical  evaluation  of  the 
detectors.  The  light  pipe,  which  has  a  total  length  of  about  3  feet,  is  made  of  |-inch  diameter 

gold  plated  thin  wall  stainless  steel  tubing. 


INCIDENT  RADIATION 
FROM  SPECTROMETER 
O.R  BLACKBODY 


WARM  BLACK 
POLYETHYLENE  FILTER 
(and  other  warm  filters 
os  required) 


Au- PLATED 
SS  LIGHT  PIPE 


COLD 

FILTER 


INTEGRATING 

CHAMBER 


TURRET  TO  HOLD 
VARIOUS  FILTERS 


COLD 

CRYSTAL  QUARTZ 
FILTER 

GaAs  DETECTOR 

CONNECTOR 
FOR  ELECTRICAL 
LEADS 


Fig.  1-1.  Detector  mount  and  filter  arrangement  for 
spectral  and  blackbody  calibration  measurements. 


The  cone  and  spherical  integrating  cavity  for 
the  detector  are  gold  plated  electroformed  cop¬ 
per.  The  cone  is  about  5  inches  long  and  the 
integrating  sphere  is  |  inch  in  diameter  with 
a  l -inch  aperture.  The  integrating  chamber 
was  used  to  minimize  differences  in  sample 
absorption.  A  rotating  turret  positioned  be¬ 
tween  the  bottom  end  of  the  light  pipe  and  the 
entrance  to  the  detector  cone  was  used  to  insert 
cold  filters  into  the  light  path  or  to  block  off 
the  light  pipe  for  dark  resistance  measurements. 
The  light  pipe,  copper  cone  and  filter  holder 
were  not  evacuated,  and  the  detectors  and  cold 
filters  were  immersed  in  liquid  He.  A  warm 
black  polyethylene  filter  was  included  at  the 
top  of  the  light  pipe,  and  a  cold  crystal  quartz 
filter  wedged  to  about  4°  was  included  at  the 
entrance  to  the  cone.  For  the  detector  calibra¬ 
tion  measurements  a  long  wavelength  pass  filter 
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was  also  used.  This  filter  has  a  sharp  room  temperature  cutoff  at  about  150p,  and  its  trans¬ 
mittance  does  not  change  significantly  on  cooling  from  room  temperature  to  4.2°K.  This  filter 
was  used  for  the  blackbody  calibration  measurements  either  at  room  temperature  (normal  back¬ 
ground)  or  at  the  detector  temperature  (reduced  background). 

Electrical  contact  to  the  sample  was  made  by  alloying  pure  tin  to  the  epitaxial  layer  at  about 
350°C  in  a  reducing  atmosphere.  The  detector  circuit  used  for  these  measurements  consists 
of  a  cooled  load  resistor  at  the  temperature  of  the  He  bath  (Rj  =  20  megohms  at  4.2°K,  26  meg¬ 
ohms  at  1.93°K)  in  series  with  the  detector  element.  One  side  of  the  detector  was  grounded, 
and  a  preamplifier  was  connected  across  the  detector.  The  signal  and  noise  measurements  were 
made  using  a  PAR  model  112  preamplifier  and  either  a  Rrower  model  132  lock-in  signal/noise 
meter  or  a  GR-1900A  wave  analyser. 

The  4.2°K  and  1.93°K  V-I  curves  for  Sample  2  under  normal  and  reduced  background  con¬ 
ditions  are  shown  in  Fig.  1-2.  Similar  behavior  was  observed  for  the  other  samples.  The  gen¬ 
eral  features  for  all  four  conditions  are  approximately  the  same.  At  low  currents  and  voltages 
the  curves  are  linear  indicating  ohmic  behavior,  but  at  higher  currents  the  voltage  does  not 

increase  as  rapidly  as  the  current.  This  decrease  in  sample  resistance  results  mainly  from 

\ 

impact  ionization  of  the  neutral  shallow  donor  levels.  At  a  critical  field,  which  depends  on  the 

electron  mobility  and  the  compensation  ratio,  the  current  increases  several  orders  of  magnitude 

while  the  field  remains  relatively  constant.  The  breakdown  field  for  Sample  2  was  about  1.5  V/cm 

at  4 . 2 ° K .  In  samples  which  are  more  heavily  compensated,  the  breakdown  field  is  larger  than 

2 

the  sustaining  field  so  that  a  negative  resistance  region  appears.  As  shown  in  Fig.  1-2,  the  back¬ 
ground  radiation  has  a  negligible  effect  on  the  breakdown  voltage,  but  there  is  a  slight  increase 
in  the  breakdown  voltage  as  the  temperature  is  decreased  from  4.2  to  1.93°K.  This  is  expected 
because  of  the  corresponding  decrease  in  mobility. 

The  resistance  of  this  detector  in  the  ohmic  region  at  4.2°K  is  about  600  kilohms  for  normal 
background  conditions,  1.4  megohms  for  reduced  background,  and  8.6  megohms  in  the  dark.  At 


Fig.  1-2.  V-I  characteristics  for  buinple  2. 
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1.93°K  the  corresponding  ohmic  resistances  are  6.5,  28.5,  and  330  megohms.  For  the  normal 
and  reduced  background  conditions  at  both  temperatures,  the  sample  resistance  is  limited  by 
the  generation  of  free  carriers  due  to  the  background  radiation. 

The  response  of  GaAs  far  infrared  detectors  at  wavelengths  longer  than  ~214  \±  is  due  prima¬ 
rily  to  a  photothermal  excitation  process^  in  which  electrons  are  optically  excited  to  the  higher 
bound  states  of  the  donor  atoms  and  then  thermally  transferred  into  the  conduction  band.  Conse¬ 
quently,  although  the  optical  absorption  by  the  impurity  centers  does  not  change  significantly  as 
the  temperature  is  reduced  from  4.2  to  1.93°K,  the  photo-conductive  response  of  the  detectors 
at  long  wavelengths  decreases  relative  to  that  for  wavelengths  shorter  than  the  photoionization 
energy.  This  is  because  the  probability  of  thermal  transfer  to  the  conduction  band  of  an  electron 
in  an  excited  state  is  much  smaller  at  the  lower  temperatures.  This  accounts  for  the  factor  of 
~11  increase  in  resistance  of  the  detector  between  4.2  and  1.93°K,  due  mainly  to  the  smaller 
number  of  free  carriers  generated  by  the  background  radiation.  (The  mobility  for  this  sample 
decreases  by  about  a  factor  of  2  over  the  same  temperature  range.)  The  larger  relative  increase 
in  resistance  at  1.93°K  between  normal  and  reduced  background  conditions  compared  to  that  at 
4.2°K  results  because  the  filter  used  for  reduced  background  conditions  blocks  a  larger  percent¬ 
age  of  the  effective  carrier-producing  background  radiation  at  the  lower  temperature. 

The  variations  of  the  signal  and  noise  voltages  and  the  signal- to-noise  ratio  with  sample 
current  are  shown  in  Fig.  1-3.  The  signal  increases  linearly  with  current  for  low  bias  values, 
but  then  bends  over  and  decreases  for  higher  currents.  The  decrease  in  signal  is  related  to  the 
onset  of  impact  ionization  of  the  neutral  shallow  donor  levels  and  the  resulting  decrease  in  de¬ 
tector  resistance.  Similar  variation  is  observed  for  the  different  temperature  and  background 
conditions  at  currents  which  produce  approximately  the  same  sample  voltage.  The  other  sam¬ 
ples  have  nearly  the  same  signal  characteristics. 

The  noise  voltage  at  low  currents  is  nearly  constant,  and  results  from  a  combination  of  the 
preamplifier  noise  with  the  Johnson  noise  from  the  sample  and  load  resistor  in  parallel.  At 


Fig.  1-3.  Signal,  noise,  and  signal/noise  dependence  on  bias  current  for  Sample  2 
at  4.  2°K  with  reduced  background  conditions. 


3 


Section  I 


higher  currents  the  noise  increases  and  reaches  a  maximum  at  about  the  same  current  which 
produces  the  maximum  signal.  The  decrease  in  noise  at  the  higher  currents  for  this  sample  is 
about  equal  to  the  decrease  in  signal  so  that  the  signal-to-noise  ratio  is  nearly  constant  over 
more  than  a  factor  of  10  range  in  bias  current.  This  property  is  useful  because  it  allows  oper¬ 
ation  without  loss  in  sensitivity  at  the  higher  bias  currents  where  the  dynamic  resistance  of  the 
sample  is  less  than  at  low  bias.  This  is  an  advantage  for  those  applications  where  response 
speed  is  critical  since  the  dynamic  RC  time  constant  for  the  detector  is  reduced.  If  the  current 
is  increased  to  still  larger  values,  the  noise  increases  several  orders  of  magnitude  and  in  this 
range  is  apparently  associated  with  the  breakdown  mechanism. 


FREQUENCY  (cm’1) 


Fig.  1-4.  Relative  spectral  response  of  photo¬ 
conductivity  for  Sample  2  at  4.  2°K. 


The  relative  spectral  response  of  Sample  2  at  4.2°K  is  shown  in  Fig.  1-4,  and  the  response 

of  the  other  samples  showed  the  same  general  shape.  The  photoconductivity  spectral  response 

measurements  were  made  using  a  Michelson  interferometer.  An  Eppley  Golay  cell  was  used 

as  a  reference  detector  to  correct  for  the  beam  splitter  and  source  response.  The  dominant 

3 

features  of  the  spectral  response  have  been  discussed  previously. 

The  peak  responsivities  of  the  detectors  were  determined  using  a  500°K  blackbody  source 
chopped  at  260  Hz  and  the  warm  or  cold  filter  described  above.  A  simple  computer  program 
was  used  to  perform  the  required  integrations  of  the  product  of  the  blackbody  emission,  the  fil¬ 
ter  transmittance,  and  the  relative  spectral  response  of  the  detectors.  Corrections  to  the  power 
reaching  the  detectors  were  made  for  the  emission  from  the  300°K  chopper  blade,  which  had  a 
measured  emissivity  of  0.955  in  the  wavelength  region  of  interest.  The  blackbody  source  and 
optical  path  were  purged  with  dry  nitrogen  to  reduce  the  water  vapor  absorption.  The  character¬ 
istics  of  the  three  GaAs  detectors  at  4.2°K  for  the  two  different  background  conditions  are  given 
in  Table  1-2. 
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TABLE  1-2 

CHARACTERISTICS  OF  GaAs  DETECTORS 

Sample 

No. 

Operating 
Conditions 
(T  =  4.  2°K) 

Peak 

Responsivity 

l!282t,(V/W) 

NEP282H(W-Hz'1/2) 

(f  =  260  Hz) 

Detector 

Dimensions 

L(mm)  X  W(mm)  X  t(p) 

Normal 

Background 

4.  OX  104 

-12 

2. 0  X  10 

1 

Reduced 

Background 

1.4X  105 

9.0X  10"13 

4.  3  X  4. 3  X  295 

Normal 

Background 

4. 0  X  104 

-12 

1 . 3  X  10 

2 

Reduced 

Background 

2.4X  105 

-13 

5.  3  X  10  J 

6.  3  X  5.  5  X  73 

Normal 

Background 

3.  OX  104 

-12 

2.8X  10 

3 

Reduced 

Background 

9.9X  104 

-12 

1. 1  X  10 

3.  8 X  3. OX  57 

G.  E.  Stillman 
C.  M.  Wolfe 
J .  O.  Dimmock 


B.  LONG  WAVELENGTH  METAL-SEMICONDUCTOR  BARRIER  LASER 

ON  P-pbo.7i5Sno.285Te 

Laser  emission  from  evaporated  metal  barriers  on  p-Pb_  _.cSn_.  .ocTe  has  been  observed 

U./1d  0  U.Zod 

at  X  =  28  p.  The  barriers  were  fabricated  by  evaporating  a  10,000-A  thick  film  on  a  freshly  etched 
<100>  surface.  An  ohmic  back  contact  was  obtained  by  evaporating  gold.  The  laser  emission 
observed  in  forward  bias  showed  Fabry-Perot  mode  structure  extending  from  X  -  27.5  p  to  28.9  p 
indicating  an  energy  gap  of  about  =  0.043  eV.  Electron  microprobe  measurements  of  the  com¬ 
position  were  performed  on  the  laser  after  etching  away  the  evaporated  Pb  barrier.  The  com¬ 
position  PbQ  ^^^Sn^  285Te  determ*-necl  in  this  manner  and  the  energy  gap  of  E^  =  0.043  eV  are  in 

satisfactory  agreement  with  the  linearly  extrapolated  results  for  E  vs  x  of  Pb.  Sn  Te  alloys 

g  l  —x  x 

of  lower  mole  fraction  SnTe.  These  devices  were  fabricated  from  large  Bridgman-grown  crys¬ 
tals.  The  material  growth  and  annealing  results  are  discussed  in  the  following  section. 

K.  W.  Nill 
A.  R.  Calawa 
T.  C.  Harman 

C.  PREPARATION  OF  MATERIAL  FOR  METAL-SEMICONDUCTOR  BARRIER  LASERS 

Experimental  results^  for  some  isothermally  annealed  Bridgman-grown  crystals  of 
Pb1  xSnxTe  for  x  =  0.13  to  x  =  0.27  were  reported  some  time  ago.  These  annealing  experi¬ 
ments  have  now  been  extended  to  x  =*  0.3  for  Bridgman-grown  crystals  equilibrated  with  a 
metal-saturated  alloy  (Pb^  3)0  5iTeo  49*  T he  carrier  concentration  results  are  displayed 
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Fig.  1-5.  Hole  carrier  concentration  determined  at  77° K  as  a  function  of  isothermal 
annealing  temperature  for  metal-saturated  Pb^  ^Sn^Te  for  x  ^  0.  3. 

graphically  in  Fig.  1-5.  It  is  seen  that  the  hole  carrier  concentration  varies  from  2.5  x  lO^/cm^ 
for  equilibrium  at  750  °C  to  1.5  X  10^/cm^  at  400  °C.  Sample  thicknesses  were  0.025  cm.  For 
the  specimens  annealed  at  lower  temperatures,  several  annealing  temperatures  were  used  and 
total  annealing  times  of  approximately  90  days  were  employed  for  each  crystal.  The  material 
used  for  the  metal-semiconductor  long  wavelength  laser  described  in  the  previous  section  is  in¬ 
dicated  in  Fig.  1-5.  ^  ^  „ 

T.  C.  Harman 

A.  E.  Paladino 


D.  TYPE  CONVERSION  AND  n-p  JUNCTION  FORMATION  IN  Hg*  Cd  Te 
PRODUCED  BY  PHOTON  BOMBARDMENT 

We  have  found  that  proton  bombardment  can  be  used  to  create  n-type  layers  in  p-type 
Hg^  xCdxTe,  and  have  used  this  technique  to  fabricate  n-p  junction  photovoltaic  detectors.  We 
feel  that  this  technique  is  exciting,  as  it  is  difficult  to  make  junction  photovoltaic  detectors  in 
this  material  by  diffusion  techniques. 

Initial  experiments  were  done  using  a  crystal  of  Hgn  .Cd  Je,  which  was  p-type  (p  ~  2  X 
10  /cm  at  room  temperature).  Using  photoresist  as  a  bombardment  mask,  15 -mil  square  areas 
in  the  sample  were  bombarded  with  100  keV  protons.  Following  bombardment,  DC  thermal  probe 
tests  and  conductance  measurements  were  made  on  both  the  bombarded  and  unbombarded  areas. 
These  measurements  indicated  that  the  bombarded  areas  had  been  converted  to  n-type,  whereas 
the  unbombarded  areas  remained  p-type.  Ohmic  behavior  was  observed  between  contacts  on  a 
bombarded  15-mil  square,  while  diode  characteristics  were  observed  between  the  bombarded 
squares  and  the  substrate.  At  room  temperature,  zero-bias  diode  resistances  as  high  as  90kil- 
ohms  were  observed. 
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The  photovoltaic  response  of  these  diodes  has  a  peak  at  1.8  jx,  with  a  quantum  efficiency  at 
the  peak  of  ~15  percent.  The  room  temperature  detectivity  at  the  peak  of  2  x  10*°cm  \  Hz/W  is 
limited  by  the  thermal  noise  of  the  detector. 

A.  G.  Foyt 
T.  C.  Harman 
J.  P.  Donnelly 

E.  A  VIDICON  WITH  A  GaAs  SCHOTTKY  BARRIER  DIODE  ARRAY  TARGET 

A  vidicon  type  television  camera  tube  has  been  constructed  with  an  image  sensing  target 
consisting  of  an  array  of  reverse  biased  Schottky  barrier  photodiodes  on  GaAs.  The  diodes  are 

accessed  by  a  low  energy  scanning  electron  beam  as  in  a  conventional  vidicon. 

14  -3 

The  n-type  GaAs  substrate  with  net  electron  concentration  of  1  x  10  cm  was  chemically 
polished  to  a  thickness  of  approximately  7  jjl  and  the  Schottky  barrier  diodes  were  fabricated  by 
electroplating  gold  or  platinum  through  photolithographically  defined  holes  in  a  pyrolytic  silicon 
oxy-nitride  layer.  The  diodes  are  reverse  biased  sufficiently  to  totally  deplete  thin  GaAs,  and 
hence  lateral  diffusion  of  the  photogenerated  carriers  is  not  a  factor  in  determining  the  resolution 
of  the  tube.  The  spatial  resolution  is  set  by  the  17.5  jjl  center-to-center  spacing  of  the  diodes 
and  the  diameter  of  the  scanning  electron  beam.  The  plated  Schottky  barriers  are  designed  to 
shield  the  passivating  silicon  oxy-nitride  layer  and  in  this  way  the  need  for  a  " resistive  sea"  is 
eliminated.  Typical  dark  currents  for  the  array  are  less  than  50  nA. 

The  GaAs  Schottky  barrier  diode  array  vidicon  exhibits  freedom  from  image  persistance  due 
to  photoconductive  lag,  is  not  damaged  by  exposure  to  intense  light  and  can  be  vacuum  baked  at 
temperatures  considerably  in  excess  of  200°C.‘  It  is  expected  that  the  target  will  have  a  broad 
spectral  response  and  an  effective  quantum  efficiency  greater  than  50  percent.  In  these  respects 
it  should  be  superior  to  a  conventional  photoconductive  vidicon.^^  Measurement  of  these  param¬ 
eters  is  now  in  progress.  Although  the  device  is  still  in  the  research  stage,  the  simplicity  of 
fabrication,  lack  of  need  for  a  resistive  sea  and  possible  benefits  in  long  term  performance  and 
stability  make  the  GaAs  Schottky  barrier  diode  vidicon  an  attractive  alternative  to  the  silicon 

array  tube7 

C.  E.  Hurwitz 
W.  T.  Bindley 
C.  M.  Wolfe 

F.  A  NEW  TECHNIQUE  FOR  THE  FABRICATION  OF  GaAs  SCHOTTKY 
BARRIER  FIELD  EFFECT  TRANSISTORS 

A  technique  has  been  developed  for  the  fabrication  of  GaAs  Schottky  barrier  field  effect 
transistors  in  which  the  source  and  drain  ohmic  contacts  are  on  one  side  of  a  thin  GaAs  layer 
and  the  gate  is  on  the  other.  This  technique  is  particularly  useful  for  the  fabrication  of  high 
frequency  devices  since  it  allows  a  reduction  of  the  source-drain  spacing  to  values  which  are 
difficult  or  impossible  to  achieve  using  a  conventional  single-sided  technique  in  which  the  gate 
must  be  placed  on  the  same  surface  between  the  source  and  drain.  Since  this  spacing  is  a  crit¬ 
ical  parameter  in  determining  the  high  frequency  performance  of  the  transistor,  substantial  im¬ 
provement  should  result. 

*  The  platinum  on  GaAs  Schottky  barriers  will  withstand  400°C  without  degradation.  Presently,  the  baking  temperature  is  limited  to 
about  200° C  by  the  epoxy  used  to  attach  the  GaAs  to  the  glass  or  sapphire  substrate  and  the  substrate  to  the  tube  envelope. 
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Fig.  1-7.  Experimental  and  theoretical  frequency  dependence  of  acoustic  gain 
coefficient.  Absolute  error  is  indicated  by  heavy  error  bars  at  2.43  GHz>  other 
error  bars  show  relative  error. 
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The  fabrication  of  the  transistors  begins  with  a  slice  of  bulk  GaAs,  one  surface  of  which  is 
patterned  and  etched  to  form  an  array  of  shallow  mesas.  (Each  mesa  will  eventually  form  an 
individual  island  of  GaAs  which  will  become  a  single  FET.)  This  surface  is  covered  with  pyro¬ 
lytic  silicon  dioxide,  and  openings  are  made  for  the  microalloyed  gold-tin  source  and  drain  ohmic 
contacts,  with  the  source-drain  gap  on  top  of  the  mesa.  Additional  gold  is  electroplated  over  the 
ohmic  contacts  and  extended  out  over  the  etched  region  to  provide  contacts  to  the  source  and 
drain.  The  wafer  is  mounted  with  epoxy  onto  a  sapphire  substrate  with  the  source-drain  side 
down  and  then  lapped  and  chemically  polished  until  only  the  individual  GaAs  islands,  the  source 
and  drain  contact  pads,  and  the  surrounding  silicon  dioxide  remain.  Finally,  the  gate  pattern 
is  defined  and  gold  is  electroplated  onto  the  GaAs  island  to  form  the  gate.  A  cross  sectional 
diagram  of  the  completed  device  is  shown  in  Fig.  1-6. 

6  —  3 

The  initial  results  have  been  obtained  on  transistors  fabricated  using  2X10  citi  n-type 
GaAs.  The  final  GaAs  layer  was  0.5|jt  thick;  the  source  drain  spacing,  5p;  the  gate  length,  5  p; 
and  the  channel  width  150p.  These  transistors  have  shown  normal  FET  characteristics,  with 
a  transconductance  of  5  to  10  millimhos. 

Several  advantages  accrue  from  this  new  fabrication  technique.  First,  bulk  GaAs  can  be 
used,  thus  obviating  the  need  for  difficult-to-make  thin  epitaxial  layers.  Second,  the  ability  to 
carry  out  fabrication  steps  on  both  sides  of  the  thin  semiconductor  layer  greatly  reduces  problems 
in  geometry  control  inherent  in  presently  used  techniques  for  FET  fabrication.  Finally,  the 
ability  to  reduce  source  to  drain  spacing  and  gate  capacitance  should  result  in  significant  im¬ 
provement  of  the  microwave  performance  of  GaAs  Schottky  barrier  field  effect  transistors. 

W.  T.  Lindley 

C.  E.  Hurwitz 
A.  G.  Foyt 

G.  FREQUENCY  DEPENDENCE  OF  ACOUSTOELECTRIC  AMPLIFICATION 
IN  n-GaAs 

Acoustoelectric  amplification  coefficients  have  been  measured  over  a  range  of  frequencies 

14  -3 

near  the  frequency  of  maximum  gain  in  oxygen-compensated  n-type  GaAs  (n  ~  5  x  10  cm  )  at 

g 

room  temperature  using  the  Brillouin  scattering  technique  described  previously.  A  meaningful 
and  accurate  gain  coefficient  can  be  obtained  from  this  material  although  it  is  relatively  inhomo¬ 
geneous  (the  most  uniform  samples  exhibited  50  percent  variations  in  resistivity  over  distances 
of  15  mm)  by  measuring  the  change  in  gain  per  unit  length  with  respect  to  a  change  in  applied 

field  dtv/dE.  This  quantity  is  independent  of  the  electrical  properties  of  the  material  for  low 

9 

fields  at  the  frequency  of  maximum  gain  in  the  qi  «  1  limit,  where  q  is  the  acoustic  wave-vector 
and  f  the  electronic  mean  free  path. 

The  measured  and  calculated  frequency  dependence  of  do/dE  for  [110]  shear  waves  is  shown 

1  0 

in  Fig.  1-7.  Relatively  good  agreement  with  the  exact  theory  is  obtained  over  the  accessible 

frequency  range  from  1.23  to  3.75  GHz.  Both  the  frequency  of  maximum  gain  and  the  magnitude 

of  the  gain  coefficient  were  found  to  be  about  15  percent  larger  than  given  by  the  qf  '<  1  approx - 

imation,  which  is  expected  since  at  4  GHz  qf  -  0.6.  The  experimental  accuracy  ±5  percent  is 

somewhat  better  than  the  accuracy  of  the  calculation  ±10  percent,  which  is  limited  primarily  by 

1 1 

the  uncertainty  in  the  best  reported  value  of  the  piezoelectric  constant. 

D.  L.  Spears 
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E.  MATERIALS  RESEARCH 


A.  USE  OF  HEAT  PIPE  IN  CRYSTAL  PULLING 


In  the  Czochralski  method  of  crystal  growth  from  the  melt,  a  seed  crystal  attached  to  a 
metal  or  ceramic  rod  grows  longer  as  it  is  gradually  pulled  from  the  melt  by  raising  the  seed 
rod.  In  most  cases  the  rod  is  simultaneously  rotated  in  order  to  give  the  crystal  rotational 
symmetry  and  to  remove  impurities  from  the  interface  region  by  stirring  the  melt.  It  is  advan¬ 
tageous  and  sometimes  necessary  to  cool  the  rod,  in  order  to  increase  the  temperature  gradient 
at  the  interface.  This  stabilizes  the  position  of  the  interface  and  reduces  the  likelihood  of  con¬ 
stitutional  supercooling.  Cooling  is  generally  accomplished  by  circulating  water  through  the  rod, 
but  this  requires  a  controlled  flow  of  water  and  the  use  of  rotating  seals  which  tend  to  leak. 

To  eliminate  the  disadvantages  of  water  cooling,  we  have  constructed  a  seed  rod  for  our 
1  2  3 

tri-arc  crystal  puller  which  is  cooled  by  means  of  a  heat  pipe,  *  a  device  with  extremely  high 
effective  thermal  conductivity.  The  heat  pipe  is  a  hollow 
evacuated  tube  containing  a  volatile  liquid  and  lined  with 
a  material  which  can  act  as  a  wick.  When  one  end  of  the 
pipe  is  heated,  the  liquid  partially  vaporizes,  the  vapor 
is  transported  at  sonic  velocity  to  the  cold  end  where  it 
condenses,  and  the  condensed  liquid  returns  to  the  hot 
end  (independently  of  gravity)  via  the  wick.  Because  the 
heat  of  vaporization  is  absorbed  at  the  hot  end  and  lib¬ 
erated  at  the  cold  end,  a  large  quantity  of  heat  can  be 
transferred  with  a  temperature  difference  between  the 
ends  which  is  typically  only  a  fraction  of  a  degree. 

The  new  seed  rod  is  shown  schematically  in  Fig.  H-l, 
in  position  for  crystal  growth.  The  copper  heat  pipe,  a 
commercial  unit  1.25  cm  in  diameter  and  30  cm  long,  has 
a  maximum  operating  temperature  of  240°C.  The  heat 
liberated  at  the  cold  end  outside  the  furnace  is  trans¬ 
ferred  to  the  room  by  means  of  a  finned  heat  exchanger 
4  cm  in  diameter  and  5  cm  long.  Because  of  the  rela¬ 
tively  high  operating  temperature,  a  teflon  seal  is  used 
at  the  furnace  gland,  as  shown  in  Fig.  II- 1 . 

The  new  seed  rod  has  been  used  to  grow  a  number  of 
single  crystals  of  VO  and  polycrystalline  boules  of  Ti  and 
Ni.  Seed  crystals  were  not  used.  Instead,  a  molybdenum 
stud  screwed  into  the  lower  end  of  the  heat  pipe  (Fig.II-1) 
was  dipped  into  the  melt  to  initiate  growth.  During  one 
run  with  Ti,  the  temperature  at  the  upper  end  of  the  heat 
pipe  was  monitored  by  means  of  a  thermocouple.  The 


Fig.  11-1 .  Seed  rod,  cooled  by  means  of 
heat  pipe,  for  tri-arc  crystal  puller. 
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temperature  was  210°  to  235°C  when  the  stud  was  in  direct  contact  with  the  superheated  melt 
(~  1 700 °C )  and  dropped  continuously  during  growth,  reaching  120°C  when  the  Ti  boule  was  about 
25  mm  long  and  8  mm  in  diameter. 

In  a  calibration  experiment  the  heat  transfer  properties  of  the  heat  pipe  were  measured  in 
still  and  moving  air.  The  lower  5  cm  of  the  pipe  was  inserted  into  a  soldering  iron  in  place  of 
the  soldering  tip,  and  asbestos  insulation  was  wrapped  around  the  iron.  The  temperature  of  the 
pipe  just  below  the  heat  exchanger  was  measured  with  a  chromel-alumel  thermocouple  wrapped 
around  the  pipe  and  covered  with  asbestos.  Figure  II- 2  shows  this  temperature  as  a  function  of 


Fig.  11-2.  Temperature  of  heat  pipe  as  a  function  of  input  power,  for  heat  exchanger 
in  still  and  moving  air. 

input  power  to  the  soldering  iron  for  still  air  and  for  a  stream  of  air  produced  by  a  14-W  muffin 
fan  about  10  cm  from  the  pipe.  In  still  air,  the  power  dissipation  is  120W  at  the  maximum  oper¬ 
ating  temperature  of  240 °C.  For  moving  air,  the  highest  available  input  power  (250  W)  gave  a 
temperature  of  only  about  150°C,  and  extrapolation  of  the  experimental  curve  indicates  that  at 
240 °C  the  power  dissipated  would  be  about  600  W,  probably  more  than  required  for  any  normal 
crystal  growth  experiment.  At  a  power  input  of  250  W,  the  temperature  difference  along  10  cm 
of  the  heat  pipe  was  only  9°C.  Under  the  same  conditions,  the  temperature  difference  along  a 
copper  cylinder  of  the  same  diameter  would  be  480°C. 

T.  B.  Reed 

B.  VAPOR  GROWTH  OF  CdO  SINGLE  CRYSTALS 

All  samples  of  CdO,  the  only  II- VI  compound  with  rocksalt  structure,  are  n-type.  Generally 

19  20  -  3 

they  have  carrier  concentrations  in  the  range  10  —  1 0  cm  ,  which  are  due  to  doubly  ionized 
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native  donor  defects  associated  with  the  presence  of  Cd  in  excess  of  the  stoichiometric  composi¬ 
tion.4  The  carrier  concentration  can  be  reduced  to  ~10*8cm  ^  by  annealing  at  elevated  temper- 

14  -3 

atures  in  an  oxygen  atmosphere,  and  values  as  low  as  ~10  cm  have  been  reported  for 

5 

sputtered  films  doped  with  copper. 

Experimental  investigations  of  CdO  have  been  hindered  by  the  high  CdO  carrier  concentra¬ 
tion  and  also  because  its  high  melting  point  (probably  above  1600°C)  and  reactivity  with  other 
oxides,  especially  silica,  make  single  crystal  growth  difficult.  Consequently  most  studies  of 
optical  and  electrical  properties  have  been  made  either  on  sintered  compacts  or  on  sputtered 
thin  films.  Crystals  approaching  1  cm  in  their  largest  dimension  have  been  grown  from  the 
vapor  phase  by  sublimation  of  CdO  powder  either  in  an  open  system  utilizing  an  inert  carrier 
gas^  or  in  a  closed  system.7  However,  these  methods  require  temperatures  exceeding  1150°C 
for  appreciable  growth  rates,  and  at  these  temperatures  the  reaction  with  container  materials 
is  severe.  Crystals  have  also  been  grown  from  CdO  powder  by  iodine  vapor  transport  at  only 

A 

about  900  to  950 °C,  but  this  method  uses  a  sealed  ampoule  made  of  fused  silica,  which  reacts 
strongly  with  CdO  even  at  this  temperature.  In  addition,  iodine  would  be  expected  to  act  as  a 
donor  impurity  if  it  were  incorporated  into  the  CdO  lattice. 

We  have  developed  an  open  tube  method  for  the  vapor  growth  of  CdO  crystals  from  the  ele¬ 
ments.  This  method  yields  sizeable  crystals  of  rather  high  purity  at  temperatures  well  below 
those  required  for  sublimation  of  CdO.  The  apparatus  is  shown  schematically  in  Fig.  II-3.  Its 


lll-l-ITItf 


principal  components  are  two  concentric  tubes  heated  by  two  adjacent  resistance  furnaces,  the 
second  hotter  than  the  first.  The  inner  tube  consists  of  a  fused  silica  section  inside  the  first 
furnace  and  a  section  made  of  high  purity  Al^O^  which  ends  at  the  center  of  the  second  furnace 
The  outer  tube  is  made  of  fused  silica,  and  the  portion  inside  the  second  furnace  contains  a  cylin¬ 
drical  liner  made  of  high  purity  Al^O^.  Only  the  two  alumina  parts  are  exposed  to  CdO  at  high 
temperature.  Attempts  to  use  alundum,  magnesia,  thoria,  and  boron  nitride  instead  of  alumina 
were  unsuccessful,  since  all  these  materials  reacted  strongly  with  CdO  at  high  temperature. 

During  a  growth  run  a  stream  of  helium  gas  (8  0  to  lOOml/min),  purified  by  being  passed 
over  titanium  metal  at  8  00°C,  enters  the  inner  tube.  It  is  then  saturated  with  cadmium  vapor 
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by  being  passed  over  a  fused  silica  boat  initially  containing  about  lOOg  of  molten  cadmium 
(nominal  purity  99.9999  percent)  heated  by  the  first  furnace  to  460°C,  where  the  vapor  pressure 
of  cadmium  is  5Torrs.  A  stream  of  pure  oxygen  (20  to  30ml/min)  enters  the  outer  tube  and 
mixes  with  the  cadmium-helium  stream  inside  the  second  furnace,  whose  center  5  inches  are 
maintained  at  1050°C  with  a  temperature  gradient  of  less  than  5°C/in.  The  CdO  crystals  grow 
downstream  from  the  end  of  the  inner  tube.  At  the  beginning  of  each  run  the  oxygen  flow  rate 
is  adjusted  while  the  growth  process  is  being  observed  through  a  window  at  the  end  of  the  outer 
tube.  Too  high  a  rate  causes  CdO  to  deposit  inside  the  inner  tube,  which  would  eventually  be 
blocked,  and  too  low  a  rate  causes  deposition  of  elemental  cadmium  at  the  downstream  end  of 
the  liner.  After  the  optimum  flow  rate  is  established,  the  run  is  continued  until  the  cadmium 
boat  is  empty,  which  requires  about  24  hours. 

Most  of  the  CdO  crystals,  which  are  black  and  quite  highly  reflective,  are  in  the  form  of 
thick  square  plates  adhering  to  the  alumina  liner  by  one  of  their  large  faces.  The  largest  crys¬ 
tals  are  formed  within  2  to  3  inches  from  the  end  of  the  inner  tube.  Several  with  dimensions 
between  5x5X2  mm  and  8x8x3  mm  are  obtained  from  most  runs.  A  number  of  intergrown 
and  twinned  crystals  with  multiple  facets  frequently  grow  from  the  tip  of  the  inner  tube,  and 
smaller  crystals  are  obtained  in  the  cooler  region  near  the  end  of  the  furnace. 

X-ray  diffraction  analysis  confirms  that  the  crystals  have  the  rocksalt  structure.  Typical 
impurity  concentrations  determined  by  mass  spectrographic  analysis  are:  Al,  about  60ppma; 

Na  and  Cl,  each  about  lOppma;  other  impurities  expected  to  be  electrically  active,  less  than 
1  ppma  (the  limit  of  detection).  The  concentration  of  Al  is  not  uniform.  Presumably  it  is  high¬ 
est  where  the  crystals  adhered  to  the  alumina  liner. 

A  few  measurements  of  the  electrical  resistivity  (p)  and  Hall  coefficient  (R^)  have  been 

made  at  room  temperature  on  as-grown  and  annealed  CdO  crystals  about  1  mm  thick.  The  car- 

20  20-3 

rier  concentrations  (n  =  —  l/R^e)  are  1.4  x  10  and  1.5  x  10  cm  for  two  as-grown  samples, 

2  0-3  ** 

1.1  x  10  cm  for  a  sample  annealed  90  hours  atl000°C  and  an  oxygen  pressure  of  about  3.8  atm, 
19  19-3 

and  2.7  x  10  and  1.8  x  10  cm  for  two  samples  from  different  growth  runs  annealed  in  air  at 

740°C  for  24  and  28  days,  respectively.  The  latter  two  values  are  significantly  higher  than  the 

19  -3 

concentration  of  1.3  x  10  cm  for  the  same  temperature  and  oxygen  pressure  interpolated  from 

4 

the  recent  rather  extensive  data  of  Koffyberg  for  the  annealing  of  CdO  crystals  grown  by  iodine 

vapor  transport.  The  Hall  mobilities  (|jl  =  R^/p)  for  the  same  five  samples  are  61,  65,  7  9, 

2  2  2  ^ 

1.4  x  10  ,  and  1.8  x  10  cm  /V-sec,  respectively.  The  latter  two  values  are  in  good  agreement 

4  19-3 

with  those  of  Koffyberg,  whose  data  do  not  extend  above  a  carrier  concentration  of  6  x  10  cm 

F.  T.  J.  Smith 

C.  STANDARD  FREE  ENERGY  CHARTS  FOR  BINARY  SELENIDES 
AND  TELLUR  IDES 

The  standard  free  energy  of  formation,  AG^.0,  is  one  of  the  most  useful  thermodynamic 
quantities  for  chemical  calculations.  In  order  to  provide  a  convenient  source  of  AG^°  data,  we 
are  preparing  an  atlas  of  charts  in  which  AG^°  is  plotted  over  a  wide  range  of  temperatures  for 
a  large  number  of  binary  compounds.  For  the  hydrides,  nitrides,  oxides,  halides,  and  sulfides, 
the  charts  have  been  compiled  from  those  scattered  through  the  literature.  For  the  selenides 
and  tellurides,  however,  no  charts  of  this  type  have  been  published,  and  it  has  been  necessary 
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to  construct  them  from  the  numerical  data  for  the  individual  compounds.  In  order  to  make  these 
new  charts  available  before  the  complete  atlas  is  published,  they  are  given  here. 

For  a  compound  MX^  (where  X  is  Se  or  Te,  M  is  any  other  element,  and  q  is  the  ratio  of 
small  whole  numbers  in  the  stoichiometric  formula),  at  a  given  temperature  AG^0  is  by  defini¬ 
tion  the  free  energy  change  accompanying  the  reaction 

-  M  +  X,  -  —  MX  (II— 1 ) 

q  2  q  q  '  ' 

when  each  reactant  is  in  its  standard  state,  i.  e.,  its  state  at  that  temperature  and  a  pressure 
of  1  atmosphere.  (The  reaction  is  written  in  terms  of  X^  rather  than  X  because  Se  and  Te  are 
diatomic  gases  over  most  of  the  temperature  range  covered  by  the  charts.)  Within  any  tempera¬ 
ture  interval  in  which  the  reactants  do  not  change  state,  AG^°  varies  linearly  with  temperature 
according  to  the  equation  AG^#  =  AH^.  —  TAS^.,  where  AH^.  and  AS^  are  respectively  the  enthalpy 
and  entropy  of  formation.  (Although  the  latter  quantities  are  weakly  temperature  dependent, 
their  average  values  within  each  temperature  interval  represent  AG^°  with  a  precision  of  better 
than  1  kcal/mole,  while  the  probable  accuracy  of  the  experimental  AG^°  data  is  about  3kcal/mole.) 
Therefore  each  plot  of  AG^°  vs  T  consists  of  a  series  of  straight-line  segments,  intersecting 
at  the  melting  and  boiling  points  of  the  various  reactants,  whose  slopes  differ  by  the  entropies 
of  fusion  or  vaporization. 

The  AG^°  charts  which  have  been  constructed  for  the  selenides  and  tellurides  are  given  in 
Figs.  II-4  and  II-5,  respectively,  and  the  values  of  AHf  and  AS^.  are  listed  in  Tables  II- 1  and  II-2, 
respectively.  The  standard  states  for  each  reactant,  whether  solid,  liquid,  or  gas,  are  indicated. 

The  quantity  AG^  is  related  to  the  equilibrium  constant  for  Eq.(II-l)  and  to  the  partial  pres¬ 
sures  p^  of  the  gaseous  species  by 


AGf  =  —  RT  In  K  =  —  HT  In 


2/q 

)2/q 


By  convention  the  pressure  is  taken  to  be  1  (unit  activity)  for  any  reactant  not  in  the  gaseous 
state.  The  quantity  AG^  is  related  to  the  electromotive  force  E  of  an  electrochemical  cell 
corresponding  to  Eq.  (II- 1 )  by  AG^  =  —  23.06  nE,  where  n  is  the  number  of  electrons  involved  in 
the  reaction,  four  in  this  case. 

The  AG^  equations  for  which  AH^  and  AS^  values  are  listed  in  Tables  II-l  and  II- 2  are  appli¬ 
cable  at  pressures  other  than  1  atmosphere.  In  each  case  the  equation  to  be  used  is  the  one  for 
which  the  standard  states  are  the  same  as  those  of  the  reactants  at  the  experimental  pressure 
and  temperature,  even  if  this  temperature  is  outside  the  range  specified  for  the  equation  in  the 

table-  T.B.  Reed 


D.  DEVIATIONS  FROM  STOICHIOMETRY  IN  EuO 

A  limited  number  of  earlier  observations  on  EuO  crystals  grown  by  directional  freezing 
from  Eu-O  solutions  of  various  compositions  indicated  that  solid  EuO  could  exist  over  a  small 

g 

but  significant  range  of  compositions.  As  the  oxygen  content  of  the  growth  solution  was  increased. 
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Fig.  11-4.  Chart  of  standard  free  energy  of  formation  (AG^°)vs  temperature  for  binary  selenides. 
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Fig.  11-5. 


Chart  of  standard  free  energy  of  formation  (AG^°)  vs  temperature  for  binary  tellurides. 
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TABLE  ll-l 

FREE  ENERGY  OF  FORMATION  OF  SELENIDES 

Compound 

Temperature 

Range 

(°K) 

-AH 

(kcal/mole) 

-AS 

(cal/°mole) 

Phase  Change  at 

Lower  Temperature* 

ZnSe 

0-  490 

99.6 

20.2 

_ 

490-  693 

102.2 

25.5 

SeM 

693-  952 

105.7 

30.6 

ZnM 

952-1184 

118.3 

43.8 

SeB 

1184-1803 

172.9 

90.2 

ZnB 

1803-2000 

151.3 

78.2 

ZnSeM 

CdSe 

0-  490 

86.3 

22.0 

— 

490-  594 

88.9 

27.3 

SeM 

594-  952 

92.0 

32.4 

CM 

952-1038 

104.6 

45.6 

SeB 

1038-1541 

152.4 

91.6 

CdB 

1541-1680 

133.9 

79.6 

CdSeM 

HgSe 

0-  490 

40.4 

21.5 

_ 

490-  630 

43.0 

26.8 

SeM 

630-  952 

71.2 

71.6 

HgB 

952-1072 

83.8 

84.8 

SeB 

1072-1200 

70.9 

72.8 

HgSeM 

GeSe 

0-  490 

39.4 

-2.6 

_ 

490-  952 

42.0 

2.7 

SeM 

952-1195 

54.6 

15.9 

SeB 

1195-1213 

9.0 

-22.2 

GeSeS  1 

1213-1500 

26.6 

-7.7 

GeM 

SnSe 

0-  490 

32.8 

-0.2 

_ 

490-  505 

35.4 

5.  1 

SeM 

505-  952 

38.8 

11.8 

SnM 

952-1133 

51.4 

25.0 

SeB 

1133-1506 

37.8 

13.0 

SnSeM 

1506-1700 

0.2 

-12.0 

SnSeB 

PbSe 

0-  490 

36.0 

2.2 

_ 

490-  600 

38.6 

7.5 

SeM 

600-  952 

40.9 

11.3 

PbM 

952-1349 

53.5 

24.5 

SeB 

1349-1531 

37.3 

12.5 

PbSeM 

1531-1700 

-0.5 

-12.2 

PbSeB 

Ag0Se 

0-  490 

20.6 

-10.7 

_ 

z 

490-  952 

23.2 

-5.4 

SeM 

952-1153 

35.8 

7.8 

SeB 

1153-1200 

21.9 

-4.2 

AgSeM 

CuSe 

0-  490 

18.9 

-9.2 

_ 

490-  952 

21.5 

-3.9 

SeM 

MnSe 

0-  490 

-56.0 

-7.9 

_ 

490-  952 

-58.6 

-2.6 

SeM 

952-1200 

-71.2 

10.6 

SeB 

*  M  —  melts 

B  —  boils 

S  —  sublimes 
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TABLE  11-2 

FREE  ENERGY  OF  FORMATION  OF  TELLURIDES 

Temperature 

Range 

-AH 

-AS 

Phase  Change  at 

Compound 

(°K) 

(kcal/male) 

(cal/0  male) 

Lower  Temperature* 

ZnTe 

0-  693 

60.9 

9.2 

_ 

693-  723 

64.4 

14.3 

ZnM 

723-1184 

72.8 

25.9 

TeM 

1184-1271 

127.4 

72.3 

ZnB 

1271-1576 

152.4 

92.0 

TeB 

1576-1700 

133.5 

80.0 

ZnTeM 

CdTe 

0-  594 

52.4 

6.9 

_ 

594-  723 

55.5 

12.0 

CdM 

723-1038 

63.9 

23.6 

TeM 

1038-1271 

111.7 

69.6 

CdB 

1271-1378 

136.7 

89.3 

TeB 

1378-1700 

120.2 

77.3 

CdTeM 

HgTe 

0-  630 

21.7 

9.3 

- 

630-  723 

49.9 

54.1 

HgB 

723-  943 

58.3 

65.7 

TeM 

943-1271 

47.0 

53.7 

HgTeM 

GeTe 

0-  723 

16.0 

1.2 

- 

723-  998 

24.4 

12.8 

TeM 

998-1213 

12.4 

0.8 

GeTeM 

1213-1215 

30.0 

15.3 

GeM 

1215-1271 

-6.0 

-14.3 

GeTeB 

1271-1700 

19.0 

5.4 

TeB 

SnTe 

0-  505 

28.4 

0 

_ 

505-  723 

31.8 

6.7 

SnM 

723-1063 

40.2 

18.3 

TeM 

1063-1271 

27.4 

6.3 

SnTeM 

1271-1530 

52.4 

26.0 

TeB 

1530-1700 

12.5 

-0.1 

SnTeB 

PbTe 

0-  600 

33.2 

2.2 

_ 

600-  723 

35.5 

6.0 

PbM 

723-1190 

43.9 

17.6 

TeM 

1190-1271 

29.6 

5.6 

PbTeM 

1271-1464 

54.6 

25.3 

TeB 

1464-1700 

15.4 

-1.5 

PbTeB 

Ag.Te 

0-  723 

14.0 

-6.4 

- 

£ 

723-1230 

22.3 

5.2 

TeM 

1230-1233 

0.2 

-12.8 

Ag2TeM 

1233-1271 

5.1 

-8.5 

AgM 

1271-1600 

30.1 

11.2 

TeB 

MnTe 

0-  723 

45.0 

-5.8 

_ 

723-1271 

53.6 

5.8 

TeM 

1271-1517 

78.6 

25.5 

TeB 

1517-1800 

85.0 

29.7 

MnM 

CoTe 

0-  723 

130.0 

5.1 

- 

723-1116 

138.4 

16.7 

TeM 

1116-1271 

142.4 

20.3 

CaM 

1271-1756 

167.4 

40.0 

TeB 

1 756-2000 

239.4 

81.0 

CuB 

sb,Te 

0-  723 

9.0 

-1.7 

_ 

£  J 

723-  893 

17.4 

9.9 

TeM 

893-1271 

8.5 

-0.1 

SbM 

1271-1600 

33.5 

24.9 

SbB 

BuTe~ 

0-  544 

12.5 

-1.3 

_ 

£  J 

544-  723 

15.9 

5.1 

BIM 

723-  860 

24.3 

16.7 

TeM 

860-1271 

15.7 

6.7 

Bi2Te3M 

1271-1600 

40.7 

26.4 

TeB 

*  M  —  melts 

B  —  boils 

S  —  sublimes 
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there  was  an  increase  in  the  oxygen  content  of  the  crystal  measured  by  combustion  analysis 
(with  x  in  EuO  increasing  from  a  minimum  value  of  0.970  for  an  Eu-rich  solution  to  a  maximum 

^  o 

value  of  1.017  for  an  O-rich  solution)  and  a  decrease  in  the  lattice  parameter  (from  5. 1450  A  to 

o 

5. 1423  A).  In  addition,  the  electrical  resistivity  of  the  crystals  varied  with  solution  composition 

and  crystal  position  in  the  ingot.  The  resistivities  of  the  first-to-freeze  portions  were  too  high 

_  \ 

to  be  measured,  while  those  of  the  last-to-freeze  portions  ranged  from  10  ohm-cm  for  crystals 
grown  from  Eu-rich  solutions  to  about  50  ohm-cm  for  crystals  grown  from  solutions  expected  to 
be  congruently  melting. 

In  order  to  obtain  further  data  on  the  influence  of  deviations  from  stoichiometry  on  the  elec¬ 
trical  properties  of  EuO,  samples  with  compositions  on  the  Eu-rich  and  O-rich  limits  of  the 
homogeneity  range  have  been  prepared  by  annealing  high-resistivity  samples  from  the  first-to- 
freeze  portions  of  a  crystal  grown  from  a  Eu-rich  solution.  (The  total  impurity  content  of  the 
crystal  detected  by  mass  spectrographic  analysis  was  less  than  200ppma.)  Each  sample  was 
annealed  at  1400°C  for  48  hours,  in  a  sealed  tungsten  capsule,  in  the  vapor  over  a  mixture  of 
two  solid  phases:  Eu-EuO  for  Eu  saturation  and  EuO-Eu^O^  for  O  saturation.  Since  in  each 
case  three  phases  (including  the  gas  phase)  are  present,  according  to  Gibbs'  phase  rule  the  par¬ 
tial  pressures  of  Eu(g)  and  02(g)  are  fixed  by  the  temperature  only.  In  order  to  achieve  equilib¬ 
rium  with  the  vapor,  the  sample  must  change  its  composition  by  solid  state  diffusion  until  it 
becomes  either  as  Eu-rich  or  as  O-rich  as  possible  at  the  annealing  temperature,  depending  on 
whether  excess  Eu  or  excess  Eu^O^  is  present  in  addition  to  EuO. 

The  annealed  samples  differ  markedly  from  the  as-grown  crystal,  whose  resistivity  was  too 
high  to  measure  at  any  temperature,  and  from  each  other.  The  resistivity  of  the  O-saturated 
sample  is  about  5  X  10^  ohm-cm  at  room  temperature  and  increases  slowly  and  smoothly  with 

4 

decreasing  temperature  to  about  5X10  ohm-cm  at  50 °K.  The  resistivity  of  the  Eu-saturated 

2  3 

sample  is  about  10  ohm-cm  at  room  temperature,  rises  to  a  maximum  of  about  10  ohm-cm 

near  80°K,  and  then  drops  by  almost  6  orders  of  magnitude  to  about  10  ^ohm-cm  at  50°K.  This 

behavior,  which  is  associated  with  the  ferromagnetic  transition  in  EuO  (Curie  point  of  69°K),  has 

also  been  observed  in  some  as-grown  samples  taken  from  the  last-to-freeze  portions  of  solution- 

q 

grown  ingots,  and  has  been  discussed  previously. 

The  striking  differences  between  the  Eu-saturated,  as-grown,  and  O-saturated  crystals 
show  that  deviations  from  stoichiometry  are  of  great  importance  in  determining  the  properties 
of  EuO,  although  no  detailed  mechanism  will  be  proposed  at  this  time.  The  earlier  results  of 
combustion  analysis  on  as-grown  samples  indicate  that  the  homogeneity  range  extends  on  both 
sides  of  the  stoichiometric  composition,  but  even  this  conclusion  should  be  regarded  as  tentative 
until  the  annealed  samples  are  analyzed.  It  should  be  noted  that  lattice  parameter  measurements 
on  these  samples  are  apparently  inconsistent  with  those  made  earlier  on  as-grown  samples.  The 
parameters  measured  for  the  Eu-saturated,  as-grown,  and  O-saturated  samples  are  respec¬ 
tively  5.1422,  5.1432,  and  5.1438A,  with  an  estimated  uncertainty  of  +0.0005A.  Thus  the  param¬ 
eter  increases  with  increasing  oxygen  content,  whereas  exactly  the  opposite  trend  was  observed 
with  as-grown  samples. 

The  method  used  to  prepare  Eu-saturated  and  O-saturated  EuO  is  based  on  the  mutual  de¬ 
pendence  of  the  compositions  of  solid  EuO  and  the  vapor  phase  with  which  it  is  in  equilibrium. 

At  any  temperature,  solid  EuO  is  stable  only  over  a  limited  range  of  Eu(g)  partial  pressures. 

The  maximum  possible  Eu(g)  pressure  is  that  in  equilibrium  with  Eu-saturated  EuO,  and  the 
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minimum  possible  pressure  is  that  in  equilib¬ 
rium  with  O- saturated  EuO.  The  available 
data  on  these  limiting  pressures  are  summa¬ 
rized  by  the  log  p-  l/T  plot  shown  in  Fig.  II-6. 
The  pressures  of  Eu(g)  are  given  by  the  upper 
loop,  whose  upper  leg  should  give  the  pressure 
over  Eu-saturated  EuO.  The  solid  line  shown 

is  actually  the  vapor  pressure  of  pure  liquid 
1 0 

Eu  metal,  which  is  given  by  log  pEu(atm)  = 
4.22  —  7,900/T.  This  represents  an  upper 
limit  for  the  pressure  over  EuO  over  the 
entire  temperature  range,  and  over  most  of 
the  range  it  should  be  a  very  good  approxima¬ 
tion  to  the  actual  Eu-saturated  pressure,  since 
the  solubility  of  oxygen  in  Eu  is  small  except 
at  temperatures  approaching  the  melting  point 
of  EuO. 

The  lower  leg  of  the  Eu(g)  pressure  loop 
in  Fig.  II- 6  gives  the  pressures  in  equilibrium 


ioVt#k 

Fig.  11-6.  Pressure-temperature  diagram  for  the  partial 
pressures  of  Eu(g)  and  02(g)  in  equilibrium  with  solid 
EuO  at  the  limits  of  its  homogeneity  range. 


withO-saturated  EuO.  The  solid  portion  gives 
1 1 

recent  data  for  these  pressures,  obtained  by  measurements  on  EuO-Eu^04  mixtures,  which 
are  represented  by  the  equation  log  Ppu(atm)  =  6.25  —  16,590/T.  The  loop  was  completed  by 
extrapolating  this  leg  and  the  upper  leg  to  meet  at  the  melting  point  of  EuO,  201 5°C  (Ref.  8). 

The  lower  loop  in  Fig.  II-6  gives  the  limiting  O 2(g)  pressures  in  equilibrium  with  EuO.  In 
this  case  the  upper  leg  is  for  O-saturated  and  the  lower  leg  for  Eu-saturated  EuO.  The  pres¬ 
sures  shown  were  calculated  from  the  Eu(g)  pressures  of  the  upper  loop  by  using  the  equilibrium 
constant  log  K  =  log  (P^Pq^)  =  18-4  ~  80,200/t,  which  was  evaluated  from  published  data  for  the 
free  energy  of  formation  of  EuO  (Ref.  12). 

It  is  seen  from  Fig.  II-6  that  at  any  fixed  temperature  the  maximum  Oz(g)  pressure  in  equi¬ 
librium  with  EuO  is  many  orders  of  magnitude  less  than  the  minimum  Eu(g)  pressure.  There¬ 
fore  EuO  has  no  congruently  subliming  composition.  If  a  sample  is  heated  in  a  dynamic  vacuum 
or  an  inert  gas  stream  it  will  lose  Eu  until  it  reaches  the  O-saturated  boundary  of  the  homogene¬ 
ity  range,  after  which  continued  vaporization  will  result  in  the  formation  of  Eu^O^. 

T.  B.  Reed 
R.  E.  Fahey 
M.  R.  Oliver 


E.  COMPOSITION  DEPENDENCE  OF  ENERGY  GAP  IN  ZnTe,  Se 
AND  ZnTe1_xSx  ALLOYS 

The  II- VI  compounds  ZnTe  and  ZnSe  form  a  complete  series  of  solid  solutions.  On  the  basis 
of  diffuse  reflectance  measurements  at  room  temperature  on  alloy  powder  samples,  Larach, 
Shrader,  and  Stocker  reported  a  minimum  in  the  energy  gap  as  a  function  of  composition  at 
about  35  mole  percent  ZnSe.  Although  later  reflectivity  measurements  by  Halsted,  Marple, 
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Fig.  11-7.  Absorption  coefficient  (a)  at  room 
temperature  vs  photon  energy  for  ZnTe  and 

ZnTe0. 674Se0. 326' 


Fig.  11-8.  Absorption  coefficient  (a)  at  77°K  vs  photon  energy  for  ZnTe^^Se^  alloys. 
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14  15 

and  Aven  and  absorption  measurements  by  Aten,  all  on  single  crystals,  supported  this  finding, 
some  workers  still  questioned  the  existence  of  a  minimum.  We  became  interested  in  this  prob¬ 
lem  because  of  our  work  with  the  oxygen  photoluminescence  band  in  these  alloys.1*"  In  order  to 
resolve  the  issue,  we  have  made  an  extensive  series  of  optical  transmission  measurements  on 
thin  samples  prepared  from  alloy  single  crystals.  Our  results  confirm  the  existence  of  the 
energy  gap  minimum. 

All  the  crystals,  including  those  of  ZnTe  and  ZnSe,  were  grown  from  the  vapor  phase  at 
temperatures  ranging  from  108  0°C  for  ZnTe  to  1180°  for  ZnSe.  Annealing  under  Zn  saturation 
at  1000°  to  1 050°C  and  furnace  cooling  generally  left  the  samples  fairly  homogeneous  and  largely 
free  of  precipitates.  Alloy  compositions  were  determined  by  electron  microprobe  analysis. 

Mass  spectrographic  analysis  of  several  boules  showed  a  few  ppma  of  a  small  number  of  impu¬ 
rities,  always  including  Cu  and  Na.  The  resistivity  of  typical  Zn-annealed  samples  ranged  from 
about  100  ohm- cm  for  ZnTe  to  10^  ohm- cm  or  higher  for  alloy  samples  containing  more  than  50 
mole  percent  ZnSe.  Samples  large  enough  to  cover  a  1/2 -mm  diameter  aperture  were  ground 
and  polished  to  thicknesses  between  20  and  7  0pm.  Transmission  measurements  were  made  with 

o 

a  grating  spectrometer  with  resolution  better  than  2 A.  Absorption  constants  in  the  region  of 

_  i 

interest,  generally  from  100  to  1000  cm  ,  were  calculated  by  using  a  value  of  0.26  for  the 
reflectivity. 

Figure  II - 7  shows  the  absorption  coefficient  (o)  at  room  temperature  as  a  function  of  photon 

energy  for  a  ZnTe  sample  and  an  alloy  sample  containing  33  mole  percent  ZnSe.  The  absorption 

_  \ 

edge  of  the  alloy,  which  has  been  measured  to  a  =  3000  cm  ,  is  seen  to  be  comparable  in  steep- 

ness  to  that  of  ZnTe.  For  both  samples  there  is  a  tail  below  about  100  cm  which  is  probably 

due  either  to  impurity  absorption,  work  damage  produced  by  polishing,  or  both.  For  ZnTe,  the 
_  \ 

energy  for  a  =  1000  cm  is  within  0.020eV  of  the  energy  gap,  which  is  known  to  be  2.25  eV  from 

1 7 

electroreflectance  measurements.  Figure  II-8  shows  similar  data  for  a  number  of  samples 
immersed  in  liquid  nitrogen.  In  general  the  absorption  edges  for  the  alloys  are  comparable  in 
steepness  to  those  for  ZnTe  and  ZnSe. 

-1 

Figure  II- 9  is  a  plot  of  the  energy  for  which  a  =  1000  cm  vs  x,  the  mole  fraction  of  ZnSe. 

At  both  77°  and  300 °K  there  is  a  minimum  in  this  energy  close  to  x  =  0.3.  The  minimum  value 
is  about  0.14  eV  below  the  value  for  ZnTe.  We  believe  that  the  energy  gaps  can  be  obtained  with 

1  3 

considerable  accuracy  by  simply  adding  0.  020  eV  to  the  plotted  values.  The  data  of  Larach,  et  al. , 
for  the  energy  gap  at  room  temperature,  obtained  by  diffuse  reflectance  measurements,  are  in 
excellent  agreement  with  our  results.  The  energy  difference  between  300°  and  77 °K  is  essentially 
constant  across  the  whole  composition  range.  It  looks  larger  near  the  minimum  but  this  is  an 
optical  illusion  due  to  the  changing  slopes  of  the  curves. 

1 8 

For  several  pseudobinary  III-V  alloy  systems  Thompson  and  Woolley  have  found  that  energy 

2 

gap  data  as  a  function  of  composition  can  be  fitted  by  the  parabolic  expression  E^,  =  A  +  Bx  +  Cx 

This  expression  also  gives  a  moderately  good  fit  to  the  present  data,  as  shown  by  Fig  II— 1 0  where 

1 9 

the  room  temperature  data  have  been  replotted.  Van  Vechten  and  Bergstresser  have  recently 

reported  that  with  one  adjustable  parameter  a  dielectric  two-band  model  in  the  virtual  crystal 

2 

approximation  can  semi- quantitatively  predict  the  coefficient  of  x  in  the  parabolic  expression 
for  the  energy  gaps  of  many  III-V  and  II-VI  alloy  systems. 
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Fig.  11-9.  Photon  energy  for  a  =  1000  cm  ^  at  77°  and  300°K  vs  mole  fraction  x  in 

ZnTe.  Se  alloys. 

1-x  x  7 


Fig.  11-10.  Photon  energy  for  a=  1000  cm  1  at  300°K  vs  mole  fraction  x  in 
ZnTe^  Se^  alloys:  comparison  of  parabolic  equation  with  experimental  data. 
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Having  confirmed  the  energy  gap  minimum  in  the  ZnTe-ZnSe  system,  we  turned  our  atten¬ 
tion  to  the  ZnTe-rich  side  of  the  ZnTe-ZnS  system,  which  is  also  of  interest  because  of  the  oxy¬ 
gen  photoluminescence  band.  The  solubility  of  ZnS  in  ZnTe  is  limited  to  about  5  mole  percent 
at  900°C  and  about  20  mole  percent  at  1200°C.20  Single  crystals  containing  up  to  15  mole  per¬ 
cent  ZnS  were  grown  from  the  vapor  phase,  although  the  yield  was  much  lower  than  for  the 
ZnTe-ZnSe  system.  The  samples  were  prepared  and  measured  in  the  same  manner  as  the 
ZnTe-ZnSe  samples  except  that  after  annealing  they  were  quenched  in  water  in  order  to  prevent 
decomposition  due  to  the  retrograde  solubility  of  ZnS. 


Fig.  11-11 .  Absorption  coefficient  (a)  at  room  temperature  vs  photon  energy 

for  ZnTe,  S  alloys. 

1-x  x  7 

Figure  II— 1 1  shows  the  room  temperature  absorption  data  for  the  ZnTe-ZnS  alloys.  Except 
for  the  sample  with  an  average  composition  of  15  mole  percent  ZnS,  the  absorption  edges  at 
sufficiently  high  absorption  coefficients  are  comparable  in  steepness  to  that  of  ZnTe.  The  15 
mole  percent  sample  was  inhomogeneous,  with  measured  composition  varying  from  about  12  to 
18  mole  percent  ZnS,  while  the  others  were  homogeneous  to  within  0.5  mole  percent. 

_  A 

In  Fig.  11-12  the  energy  for  a  =  1000  cm  is  plotted  versus  x,  the  mole  fraction  of  ZnS. 
Although  the  energy  gap  of  ZnS  is  1.5  eV  greater  than  that  of  ZnTe,  adding  ZnS  to  ZnTe  causes 
a  rather  rapid  decrease  in  the  energy  gap  of  the  alloys.  The  curve  shown  is  not  extended  to  the 
point  at  x  =  0.15  because  of  the  uncertainty  in  that  point.  A  minimum  cannot  be  observed  because 
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Fig.  11-12.  Photon  energy  for  a  =  1000  cm  at  room  temperature  vs  mole  fraction  x 

in  ZnTe,  S  alloys. 

1  -x  x  ' 

samples  with  high  enough  ZnS  content  cannot  be  prepared.  However,  it  is  clear  that  the  ZnTe- 
ZnS  system  provides  another  example  of  a  strongly  nonlinear  dependence  of  energy  gap  on  alloy 
composition.  G  w.  Iseler 

F.  ENERGY  GAP  MINIMUM  IN  CdTe1_xSex  ALLOYS 

The  composition  dependence  of  the  energy  gap  in  CdTe-CdSe  alloys  has  been  determined  by 
means  of  optical  transmission  measurements  at  room  temperature.  The  data  establish  the  exist¬ 
ence  of  an  energy  gap  minimum,  which  was  anticipated  at  the  outset  of  the  investigation  because 

1  3 

such  a  minimum  had  been  reported  in  the  ZnTe-ZnSe  system,  although  it  was  not  universally 
accepted  at  that  time.  (The  minimum  in  the  ZnTe-ZnSe  system  has  now  been  confirmed  by  the 
study  described  in  the  preceding  section  of  this  report.) 

21 

The  phase  diagram  for  the  CdTe-CdSe  system  has  recently  been  determined.  Since  CdTe 
has  the  cubic  zincblende  structure  and  CdSe  has  the  hexagonal  wurtzite  structure,  a  complete 
series  of  solid  solutions  with  the  same  structure  is  not  formed.  Instead,  there  are  two  broad 
single-phase  regions  in  the  solid  phase,  one  including  CdTe  and  alloys  with  zincblende  structure 
and  the  other  including  CdSe  and  alloys  with  wurtzite  structure.  These  are  separated  by  a  two- 
phase  region  only  about  3  mole  percent  wide,  whose  boundaries  shift  toward  increasing  CdSe 
content  with  decreasing  temperature.  For  compositions  between  about  30  and  50  mole  percent 
CdSe  it  is  possible  to  obtain  alloys  with  either  structure,  depending  upon  preparation,  annealing, 
and  cooling  conditions. 

The  samples  studied  in  this  investigation  were  single  crystals  cut  from  polvcrystalline  ingots 
grown  from  the  melt  by  directional  freezing.  Most  of  the  ingots  were  prepared  from  stoichio¬ 
metric  melts  by  the  Bridgman  method.  One  ingot  containing  about  60  mole  percent  CdSe  was 

prepared  from  stoichiometric  material  by  a  zoning  technique,  and  an  ingot  of  CdSe  was  obtained 

22 

from  an  Sc-rich  solution  by  the  temperature  gradient  solution  zoning  method.  According  to  the 
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phase  diagram,  melts  containing  more  than  20  mole  percent  CdSe  crystallize  in  hexagonal  form. 
However,  three  cubic  crystals  were  obtained  with  compositions  between  30  and  45  mole  percent 
CdSe.  One  of  these  was  transformed  to  the  cubic  form  during  furnace  cooling  of  a  Bridgman- 
grown  ingot,  and  the  others  were  transformed  from  as-grown  hexagonal  crystals  by  annealing 

i 

at  temperatures  below  the  hexagonal-to-cubic  transition  temperature.  Measurements  were  also 
made  on  a  sample  containing  49  mole  percent  CdSe  which  X-ray  diffraction  data  indicated  to  be 
a  polytype  of  unidentified  structure.  This  polytype  was  apparently  formed  during  furnace  cool¬ 
ing  as  a  metastable  intermediate  in  the  hexagonal-to-cubic  transformation. 

Transmission  measurements  in  the  vicinity  of  the  absorption  edge  were  made  at  room  tem¬ 
perature  with  a  single  pass  prism  spectrometer  on  samples  ground  to  a  thickness  of  about  8  0pm 
and  optically  polished  but  not  etched.  Absorption  coefficients  were  calculated  by  assuming  a 
fixed  reflectivity  for  each  composition,  ranging  from  0.24  for  CdTe  to  0.20  for  CdSe. 

Figure  11-13  shows  semi-logarithmic  plots  of  absorption  coefficient  vs  photon  energy  for 
representative  samples  with  cubic,  hexagonal,  and  polytype  structure.  The  compositions  of  the 
alloys,  which  were  determined  by  electron  microprobe  analysis,  are  specified  by  giving  the  mole 
fraction  of  CdSe,  that  is,  x  in  CdTe^  x^ex-  F°r  CdTe,  CdSe,  and  alloys  with  x  less  than  0.4, 
the  plots  are  linear  at  the  upper  end  of  the  absorption  coefficient  range.  For  these  alloys,  the 

absorption  edges  are  comparable  in  steepness  to  that  of  CdTe,  which  is  in  good  agreement  with 
23 

the  data  of  Marple.  The  absorption  edges  for  alloys  with  higher  CdSe  content  are  nonlinear 
and  less  steep.  In  all  cases  there  is  a  tail  at  the  lower  end  of  the  absorption  range  which  is 
probably  due  to  impurity  absorption,  work  damage  produced  during  polishing,  or  both. 

The  energy  gap  for  each  sample  was  found  by  adding  0. 020  eV  to  the  photon  energy  for  which 

-1 

the  absorption  coefficient  is  400  cm  .  This  procedure  gives  values  of  1.50  and  1.71  eV  for  the 
energy  gaps  of  CdTe  and  CdSe,  respectively,  in  good  agreement  with  the  published  values.  The 
energy  gaps  obtained  in  this  manner  for  all  the  samples  are  plotted  vs  mole  fraction  of  CdSe  in 
Fig.  H-14.  Data  have  been  obtained  for  cubic  samples  with  x  between  0  and  0.45.  Alloying 
CdTe  with  CdSe  initially  produces  a  marked  decrease  in  the  energy  gap.  This  continues  until 
about  x  =  0.35,  after  which  there  is  no  significant  change.  These  results  are  sufficient  to  estab¬ 
lish  the  existence  of  an  energy  gap  minimum  in  the  CdTe-CdSe  system.  The  density  of  data 
points  for  hexagonal  samples  is  much  lower  than  for  cubic,  because  the  preparation  of  alloys 
with  high  CdSe  content  is  made  difficult  by  their  high  melting  points.  The  limited  data  available, 
particularly  for  the  point  at  x  =  0.35,  indicate  that  the  energy  gap  is  systematically  higher  by 
about  0.03eV  for  hexagonal  samples  than  for  cubic.  A  difference  of  this  magnitude  would  not  be 
surprising,  since  the  energy  gap  of  ZnS,  which  is  about  3.8  eV  at  room  temperature,  is  about 
0.1  eV  greater  for  the  hexagonal  form  than  for  the  cubic.  The  energy  gap  value  for  the  polytype 
alloy  sample,  which  appears  to  be  intermediate  between  the  values  for  hexagonal  and  cubic  sam¬ 
ples,  is  also  consistent  with  the  results  for  ZnS,  since  the  energy  gaps  of  ZnS  polytypes  lie  be- 

24 

tween  those  of  the  wurtzite  and  zincblende  forms. 

The  energy  gap  data  for  the  cubic  and  hexagonal  CdTe-CdSe  alloys  as  a  function  of  x  are 

well  represented  by  the  curves  in  Fig.  11-14  labelled  (1)  and  (2),  respectively.  These  curves 

were  calculated  from  the  equations  shown,  which  have  the  parabolic  form  used  by  Thompson  and 
1 8 

Woolley  to  fit  energy  gap  data  for  a  number  of  III-V  alloy  systems.  The  same  coefficients  of 
2 

x  and  x  have  been  chosen  for  both  equations,  so  that  the  only  difference  between  them  is  in  the 
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Fig.  11-13.  Absorption  coefficient  (a)  at  room  temperature  vs  photon  energy 

for  CdTe.  Se  alloys. 

1  -x  x  1 
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value  of  the  first  term,  which  gives  the  energy  gap  of  CdTe.  Therefore,  the  calculated  energy 
gap  values  for  the  cubic  and  hexagonal  forms  differ  by  a  constant  amount,  0.03eV,  over  the 
entire  composition  range.  (The  apparent  variation  in  this  difference  in  Fig.  II- 14  is  an  optical 
illusion.) 

Figure  11-14  also  shows  the  energy  gap  data  for  CdTe-CdSe  alloys  obtained  by  Prytkina, 

25 

et  al. ,  from  absorption  measurements.  They  too  observe  a  minimum  in  the  energy  gap,  but 
in  comparison  with  our  results  the  minimum  is  somewhat  deeper  and  occurs  at  higher  CdSe  con¬ 
tent.  In  addition,  their  energy  gap  values  are  systematically  lower  than  ours,  as  shown  by  the 
results  for  CdTe  and  CdSe.  It  is  not  possible  to  give  an  explanation  for  these  differences,  since 
Prytkina,  et  al.,  do  not  present  any  absorption  data  or  describe  their  procedure  for  evaluating 
the  energy  gap  from  the  data.  A  j  Strauss 


G.  ENERGY-LEVEL  SCHEME  FOR  Ni^O.Li 

Evidence  for  localized  d  electrons  in  MnO,  Fe^_^0,  CoO  and  NiO  comes  from  optical- 
absorption  spectra,  the  magnitudes  of  the  atomic  moments,  and  Jahn-Teller  distortions  below 

ln-i-Hmfil 

electric  loss  measurements  on  NiO:Li  ceramics  primarily  m*  (dn4*) 


'2fl ( F«z  : d*  ) 


,  (Co* 


_  *A,fl(Mn24 


d3) 


(  N  i  .*  d  ) 


indicate  that  any  activation  energy  for  hopping 
of  the  hole  trapped  at  an  Li+  center  is  <  0.01  eV, 
and  transport  measurements  have  shown  that 
although  the  activation  energy  for  the  conductiv¬ 
ity  is  0.4  <  q  <  0.85  eV  in  the  temperature  range 
300  <  T  <  1 000°K,  the  mobilities  are  not  acti¬ 
vated.  These  facts  have  led  to  a  variety  of  spec¬ 
ulations  in  the  literature  about  the  character  of 
the  d  electrons  and  of  the  positive  charge  car¬ 
riers.  In  none  of  these  speculations  has  serious 
attention  been  given  to  the  role  of  native  defects. 

Optical  absorption  and  electroreflectance 

spectra  permit  placement  of  the  ground  state  of 

the  d- state  manifold  relative  to  the  edge  of  the 

conduction  band,  as  shown  in  Fig.  11-15.  The 
6  2+ 

A^(Mn  )  level  is  relatively  stable  because  of 
the  large  intra-atomic  exchange  stabilization 
associated  with  a  half  filled,  high  spin  d  shell. 

Otherwise  there  is  a  progressive  stabilization 
with  increasing  atomic  number.  From  band- 

calculation  estimates  and  a  shoulder  in  the  absorption  spectrum,  the  energy  gap  between  the 

2  _l  8 

O  :2p  bands  and  the  Ni  :d  4s  band  in  NiO  is  about  5.5eV. 

MnO,  Fe1  fiO,  CoO  and  NiO  are  all  anion-excess  compounds.  The  predominant  defects  are 
cation  vacancies.  Given  localized  d  electrons,  an  isolated  cation  vacancy  Vc  traps  a  pair  of 


Fig.  11-15.  Energy  level  diagram  for  MnO,  FeO, 
CoO,  and  NiO  in  the  absence  of  lattice  defects. 


29 


Section  II 

trivalent  cations,  one  on  each  side  of  it,  and  the  first  and  second  ionization  potentials  may  be 
estimated  as 

«  (2e2/«r  R)  -  (eZ/2«-R)  «  0.6eV  (R-2) 

e2  »  2e2/*R  «  0.8eV  (II-3) 

where  k  «  12  is  the  dielectric  constant  and  R  is  the  nearest-neighbor  cation-cation  separation. 
The  corresponding  modification  of  Fig.  11-15  is  the  creation  of  two  M2+:VC  levels  above  the  filled 
ground-state  IVP  level,  as  shown  in  Fig.  11-16  for  NiO.  Although  the  Ni^+:VC  levels  are  empty 
in  the  diagram,  the  axis  of  the  Ni^-V^Ni111  triplet  may  be  rotated  by  an  external  electric  field, 
the  holes  being  free  to  hop  between  the  twelve  cations  neighboring  the  cation  vacancy. 
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Fig.  11-16.  Energy  level  diagram  for  NiO, 
showing  levels  due  to  cation  vacancies  and 
Li  impurity  atoms. 


Low-spin  Ni  1  ions  are  created  from  high-spin  Ni2  ions  by  the  removal  of  an  e  electron 

2-  ° 

(electron  in  an  orbital  that  a  bonds  with  O  ions),  and  these  may  be  localized  to  the  cation  to 
create  discrete  Ni111  ions.  On  the  other  hand,  the  orbitals,  which  are  directed  toward  near¬ 
neighbor  cations,  are  strongly  perturbed  by  a  cation  vacancy.  In  fact,  each  cation  vacancy  would 
create  symmetrized  trap  orbitals  consisting  of  linear  combinations  of  t^  orbitals  from  neigh¬ 
boring  cations  admixed  with  anion  p^  orbitals  that  would  a-bond  to  the  vacancy. 

In  TiO  and  VO  ,  where  the  e  orbitals  are  empty  and  the  t~  electrons  are  itinerant,  these 
x  x  g  g 

trap  orbitals  are  lifted  above  the  Fermi  energy.  In  NiO,  where  localized  holes  are  created  by 
the  removal  of  an  e^  electron,  removal  of  a  trap  electron  requires  the  additional  crystal-field 
energy  lODq.  However,  the  added  covalent  mixing  in  the  trap  orbitals  reduces  the  effective 
1 0 Dq  by  about  1  eV,  which  places  the  energy  of  the  trap  orbitals,  designated  ^g^c  *n 

about  1  eV  below  the  Ni2+:V  localized-hole  levels,  or  approximately  at  the  3A~~(Ni2+ :d8)  level 

C  ^ 

in  Fig.  11-16.  In  addition,  a  second  trap  orbital  is  lifted  out  of  the  O  :2p  bands,  in  analogy  to 
the  anion-vacancy  trap  orbital  dropped  from  the  conduction  band  in  alkali  halides  to  form  a  color 
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center.  This  cation-vacancy  trap  orbital  is  designated  2p:V  in  Fig.  11-16  Because  the  t~  :V 

c  c 

and  2p:  trap  orbitals  have  the  same  symmetry,  they  mix  with  one  another.  This  means  that 

the  closer  the  2p:Vc  level  comes  to  the  1 2g;Vc  ^eve^  the  more  2p  character  is  imparted  to  the 

t~  -V  orbital  and  hence  the  more  extended  it  is  in  real  space.  It  follows  that,  with  strong  mix- 
£g  c 

ing  and  relatively  large  concentrations  [  V^j  of  cation  vacancies,  holes  in  the  t^V^  levels  may 
tunnel  without  activation  energy  between  neighboring  sites. 

In  Fig.  11-16,  the  t^ g‘-V  c  level  is  shown  broadened  into  a  narrow  impurity  band.  Finally,  the 
Fermi  energy  e ^  would  lie  between  the  t2g:^c  *eve^  (or  band)  and  the  first-ionization  level  at 
However,  from  Fig.  11-15  it  is  clear  that  any  iron  or  cobalt  impurities  would  act  as  donor  centers: 
Fe24  -*■  Fe^++  e  ,  where  the  electrons  fall  into  the  levels  at  Thus  the  model  accounts  well 

for  0.4  <  q  <  0.85  eV  from  the  resistivity  vs  temperature  measurements,  for  variation  in  a  vs 

2 

[P021  wi-th  Iron  impurities,  and  for  isotropic  p-type  conduction  having  a  small  »  0.3  cm  / 

V-sec),  unactivated  mobility  while  simultaneously  preserving  the  localized  character  of  the  d 
2+ 

electrons  at  the  Ni  ions  not  neighboring  a  cation  vacancy. 

2+  2—  J2J  2+  2—  2+ 

Since  the  Ni  -O  -Ni  interactions  are  ferromagnetic  whereas  Ni  -O  -Ni  interactions 

are  antiferromagnetic,  it  is  also  possible  to  rationalize  with  this  model  a  Aq  ~  0.2eV  on  raising 

the  temperature  through  T^,  an  absorption  peak  at  0.24  eV  below  T^,  and  a  photocurrent  peak 

at  0.23  eV  in  room-temperature,  amorphous  NiO  films. 

In  order  to  generalize  Fig.  II- 1 6  to  include  lightly  Li-doped  NiO,  it  is  only  necessary  to  add 

2+ 

Ni  :Li  levels,  which  represent  the  trapping  energy 


€  Li  -  e  /k  R  »  0.4  eV 


(II -4) 


of  an  Ni***  ion  to  a  near-neighbor  site  of  an  Li+  ion.  This  energy  is  in  good  agreement  with  the 
appearance  of  an  additional  optical  absorption  peak  at  0.43  eV  in  Li-doped  NiO 

J.  B.  Goodenough 


H.  POLARON  MORPHOLOGY  IN  VANADIUM  SPINELS 
FROM  SEEBECK  COEFFICIENTS 

From  a  study2^  of  the  series  of  spinels  M2+[V2+]  O^,  it  was  possible  to  infer  that  the  critical 

V-V  separation  R  for  a  transition  from  localized  to  itinerant  vanadium  d  electrons  is  approached 

c  °  +  27 

in  Mg[V2]  O^  and  Zn[V2]C>4,  where  R  «  2. 97  A.  Li  [V2]04is  metallic  and  Pauli  paramagnetic,  in¬ 
dicating  itinerant  d  electrons,  and  variations  in  the  physical  properties  with  x  in  A2_xLix[V2l  O^, 

o 

A  =  Mg  or  Zn,  have  suggested  a  critical  composition  x  «  0.55,  corresponding  to  R  «  2.94  A,  for 

c  28  29 

a  transition  from  small-polaron  to  itinerant  d  electrons.  However,  it  was  noted  that  the 
Seebeck  coefficients,  although  temperature-independent  as  required  by  the  small-polaron  formula 

e+  =  198  log  [(N-p)/p]  +  S*/e  (U-5) 

2  + 

are  not  described  by  Eq.  (II- 5 )  either  in  this  system  or  in  the  systems  Mg[Bx  w^ere 

B  =  Mg  or  Ni,  if  the  number  N  of  sites  available  to  the  p  polarons  is  taken  as  the  number  of  va¬ 
nadium  ions.  In  oxides  S*/e  <  10pV/°K  is  small  relative  to  the  first  term,  so  that  Eq.  (II- 5 )  con¬ 
tains  no  adjustable  parameters  and  therefore  represents  a  strong  constraint  on  any  microscopic 
model  for  the  polarons.  Since  p  is  given  by  the  chemical  composition,  the  model-dependent  param¬ 
eter  is  N.  Adjustment  of  N  to  allow  for  intermediate  size  polarons  has  provided  a  successful 
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TABLE  11-3 

WEIGHTED  PROBABILITY  COEFFICIENTS  an  =  DESCRIBING  1 

OF  POLARON  SITES  ASSOCIATED  WITH  A  B-SITE  VANADIUM 
HAVING  n  MAGNESIUM  NEAREST  NEIGHBORS 

NUMBER 

ION 

n 

Pn 

Clusters 

X 

n 

a 

n 

0 

1 

>* 

CM 

1 

1 

1 

6 

2V2  +  V4 

5 

2 

15 

2 

15 

|(v2+v4)+ f  <4v 

3 

45 

3 

20 

11 

4 

55 

4 

15 

2V2 

2 

30 

5 

6 

CM 

> 

1 

6 

6 

1 

V4+ 

2 

2 

Fig.  11-17. 
Mg[Mg  V9 
spinels? 


Seebeck  coefficient  vs  x  in 
x]04  and  Mg[NixV2_xJ04 


32 


Section  II 


30 

interpretation  of  the  temperature-independent  Seebeck  coefficients  in  the  systems, 

which  encourages  use  of  the  Seebeck  coefficient  to  obtain  indirect  evidence  of  the  morphology 
of  the  polarons  in  the  vanadium  spinels. 

Examination  of  the  crystal  structure  suggests  three  possible  morphologies  for  polarons  in 

spinels:  (1)  the  charge  carrier  is  localized  to  a  discrete  atomic  position,  creating  an  identifiable 
4+ 

V  ion.  (2)  The  charge  carrier  is  shared  between  two  neighboring  cations,  occupying  a  V2- 
cluster  molecular  orbital.  (3)  The  charge  carrier  is  shared  by  the  four  neighboring  cations  of 

an  elementary  rocksalt-type  cube,  occupying  a  V4-cluster  molecular  orbitaL  For  the  spinels 

2+  4 
A  [V2)04,  the  molecular  concentration  N  of  available  sites  for  each  type  of  p-type  polaron  is 

N(V4+)  =  2  ,  N(V2)=6  ,  N(V4)=1  (II-6) 

The  third  of  these  provides  the  only  possible  alternative  for  describing  with  Eq.  (U-5)  the 

2+  2+ 

measured  Seebeck  coefficients  in  A,  Li  [V~]CL,  A  =  Mg  or  Zn,  and  Mg[B  V9  ]Oa,  B  =  Mg 

or  Ni,  spinels.  However,  with  V4-cluster  polaron  sites  in  the  end  members  Mg[V~]CL  and 

4  2-f  2+  ^  4 

Zn[V2]04,  it  is  clear  that  in  the  systems  Mg[Bx  V2_x]04  each  B  ion  disrupts  the  formation 

of  V4  clusters,  forcing  nearest-neighbor  vanadium  ions  to  participate  in  clustering  instead. 

Therefore,  in  order  to  obtain  N  for  the  system,  we  must  multiply  the  probability  (1  —  of 

finding  a  vanadium  ion  on  a  particular  site  by  the  sum  of  the  probabilities  of  finding  0,  1, ...  , 

6  Mg  nearest  neighbors,  each  properly  weighted  by  a  factor  \n  that  reflects  the  polaron  clusters 

in  which  the  vanadium  ion  participates.  This  weighting  factor  counts  each  V.  cluster  as  l/2, 

4+  4 

each  cluster  as  1  and  the  discrete  V  ion  as  2.  This  leads  to 

6 

N  =  (l-f)  £  an  (f  )n  (1  -  f  )6"n  and  p  =  x  (II-7) 

n=0 

where  the  are  given  in  Table  II - 3 .  Substitution  of  Eq.  (II-7)  into  Eq.  (II- 5 )  with  the  constraint 
S*/e  =  0  gives  the  solid  curve  of  Fig.  11-17. 

Substitution  of  tetrahedral- site  Li  ions,  on  the  other  hand,  can  only  increase  the  size  of 
the  polarons  in  the  vicinity  of  the  Li+  ions,  since  the  electrons  are  all  itinerant  in  the  end  mem¬ 
ber  Li[V2]04.  The  simplest  model  would  have  large-polaron  clusters  consisting  of  all  the  vana¬ 
dium  ions  nearest-neighbor  to  Li+  ions.  A  negligible  contribution  to  the  Seebeck  coefficient 
from  charge  carriers  in  these  molecular  orbitals  is  equivalent  to  2x  polaron  sites  for  p  ~  x. 

Since  the  only  intermediate -polaron  sites  are  those  V4  clusters  having  only  A -site  Mg2+  or 
Zn2  ions  as  nearest  neighbors,  this  gives  an  "effective" 

N  =  2x  +  (1  —  x)^  and  p  =  x.  (II- 8 ) 

2  9 

Jaskowsky,  although  arguing  from  a  different  physical  model,  has  already  shown  that  these 
values  of  N  and  p  in  Eq.  (II- 5 )  give  an  excellent  description  of  the  measured  Seebeck  coeffi¬ 
cients  if  S*/e  =  0,  corresponding  to  no  adjustable  parameter. 

In  Fe[V2]04  the  V-V  separations  are  R  =  2.99A,  and  in  the  system  Fe[FexV2_x]04  a  con¬ 
ventional  model  (N  =  2  -  x)  is  applicable  for  x  >  0.15,  the  compositional  range  in  which  all  the 
iron  ions  are  divalent. 
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These  results  illustrate  use  of  a  temperature-independent  Seebeck  coefficient  to  investigate 

polaron  morphology  and  indicate  that  in  the  vanadium  spinels  the  polarons  change  with  decreas- 

4+ 

ing  V-V  separation  from  identifiable  V  cations  to  V^-cluster  and  V^-cluster  polarons,  and 

finally  to  itinerant  electrons,  in  a  narrow  range  2.99  >  R  >  2. 94 A.  Furthermore,  in  the  systems 
2+ 

A  =  Mg  or  Zn,  the  polaron  to  itinerant-electron  transition  does  not  occur  homo¬ 
geneously  at  x  «  0.55,  but  inhomogeneously  in  the  vanadium  shell  of  nearest  neighbors  to  the 
.+ .  c 

L1  lons'  J.  B.  Goodenough 


I.  CRYSTAL  STRUCTURE  OF  HEXAGONAL  RbNiF3(6H) 

31 

r  NiF3  exists  in  two  modifications.  At  ambient  pressure,  it  has  the  hexagonal  BaTi03(6H) 

struci^re-*2  while  at  high  pressure  (~25kbars)  its  structure  transforms  to  the  cubic  perovskite 

33  34  ° 

type  '  with  a  =  4. 07  A.  This  pressure  transformation  is  similar  to  those  found  for  many  com¬ 
pounds  with  a  general  formula  ABX^,  having  hexagonal  or  cubic  close-packed  AX^  layers.33  3 

The  hexagonal  form  of  RbNiF.^  is  of  considerable  interest  because  of  its  ferrimagnetic  and 
optical  properties.  Since  accurate  bond  distances  and  angles  are  necessary  for  any  quantitative 
interpretation  of  physical  properties,  we  undertook  the  refinement  of  the  atomic  position  param¬ 
eters  of  RbNiF^(6H),  as  well  as  the  determination  of  the  temperature  dependence  of  the  lattice 
parameter. 


o  o 

We  find  the  lattice  parameters  at  room  temperature  to  be  a  =  5.843 A  and  c  =  14. 309 A,  which 

31 

are  identical  to  those  previously  reported.  The  variation  of  the  lattice  parameters  in  the  tem¬ 
perature  range  -175°  to  600°C  is  shown  in  Fig.  11-18.  No  change  in  symmetry  or  anomaly  in  the 
temperature  dependence  of  the  lattice  parameters  was  noted  over  the  entire  temperature  range. 


which  includes  the  Curie  point. 


Fig.  11-18.  Crystallographic  parameters  of 
RbNiF^(6H)  as  a  function  of  temperature. 


Fig.  11-19.  Essential  features  of  the 
RbNiFg(6H)  structure. 
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RbNiF3(6H)  has  the  hexagonal  BaTiO^  structure  indicated  in  Fig.  11-19.  It  contains  face- 
shared  pairs  of  NiF^  octahedra  connected  by  single  octahedra  sharing  only  corners.  The  rubid¬ 
ium  ions,  which  fill  the  six  large  voids  per  unit  cell  of  this  octahedra  network  have  twelvefold 
coordination.  There  are  two  different  positions  for  each  type  of  atom,  as  shown  in  Fig.  11-19. 


TABLE  11-4 

CRYSTAL  STRUCTURE  OF  RbNiF^H) 

Space  Group 

P6g/mmc  (No.  194) 

Unit  Cell  Dimensions 

o  =  5. 843  A;  c=  14. 309  A 

Cell  Contents 

6  RbNiF3 

Rbj  in  (2b) 

±(0,0, 1/4) 

Rb(|  in  (4f) 

±  (1/3, 2/3,  z;  2/3, 1/3, 1/2  +  z)  z  =  0.0954(2),  0.0954(7)* 

Nij  in  (2a) 

(0,0,0;  0,0, 1/2) 

Ni(|  in  (4f) 

±  (1/3, 2/3, z;  2/3, 1/3, 1/2  +  z)  z  =  0.  8462(5),  0.  8450(8)* 

Fj  in  (6h) 

±  (x,2x,  1/4;  2x,x,3/4;  x,x,l/4)  x  =  0.517(l),  0.5023(82)* 

F,,  in  (12k) 

±(x,2x,z;  2x,x/z;  x,x,z;  x/2x,l/2  — z;  2x,x,l/2  +  z/ 

x,x,l/2  +  z)  x  =  0.  830(3),  0.8313(89)*;  z  =  0.081(l),  0.0740(17)* 

Cell  Temperature  Factor 

B.=  1.9(1) 

*  Our  results  are  listed  first  with  estimated  error  limits  in  parentheses;  Babel's  results  (37) 

are  listed  second  for  comparison. 

The  structure  has  only  five  variable  parameters.  Since  the  powder  pattern  is  not  complex, 

we  chose  to  use  intensity  data  collected  on  powdered  samples.  Refinement  based  on  thirty-five 

intensities  gave  the  cell  temperature  factor  and  the  five  variable  position  parameters  listed  in 
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Table  n-4  with  a  reliability  factor  R  =  5.8.  Babel’s  position  parameters  from  a  single -crystal 
structure  refinement  are  also  included  in  Table  E-4.  They  show  good  agreement  with  our  cation 
parameters,  but  significant  differences  with  our  fluorine  positions.  Most  affected  by  this  dis¬ 
crepancy  are  the  Fj-Fj  distances  and  Nijj-Fj-Nij ^  angle.  We  find  these  distances  to  be  2.62  and 
3.22A,  compared  to  Babel’s  values  of  2.88  and  2.96A,  and  we  find  the  angle  to  be  84.5°,  com¬ 
pared  to  Babel’s  value  of  78.5°. 

One  can  only  speculate  as  to  why  there  is  a  difference  between  our  results  and  those  of  Babel. 

We  do  feel  however  that  our  atomic  positions  are  more  probable  since  they  are  in  agreement  with 
32  37  38  39 

similar  structures,  9  *  NMR  data  on  RbNiF^(6H),  and  our  understanding  of  electrostatic  re- 

3  40 

pulsion  and  screening.  A  detailed  discussion  of  the  RbNiF^  structure  will  be  presented  elsewhere. 


R.J.  Arnott  J.  R.  O'Connor 
J.  M.  Longo 
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X-RAY 

TABLE  11-5 

DATA  FOR  ATMOSPHERIC  PRESSURE  FORM  OF  CsMnCI3(9R) 
a  =  7.  288  A,  c  =  27. 44  A  (Ref.  42) 

hki 

^calc 

*obs 

^obs 

1  0 

4 

25.5 

24 

4.64 

0  1 

5 

43.5 

45 

4.  13 

1  1 

0 

87.4 

80 

3.64 

1  0 

7 

8.6 

10 

3.33 

0  2 

1 

12.0 

13 

3.  13 

2  0 

2 

5.2 

6 

3.07 

0  0 

9 

29.4 

34 

3.04 

0  1 

8 

23.5 

20 

3.01 

0  2 

4 

71.5 

67 

2.864 

2  0 

5 

100.0 

100 

2.732 

1  0 

10 

10.3 

13 

2.516 

0  2 

7 

14.  1 

18 

2.458 

1  1 

9 

3.4 

2 

2.  342 

2  0 

8 

21.  1  1 

28 

2.319 

0  1 

11 

3.8  1 

2  1 

4 

9.5 

8 

2.252 

1  2 

5 

17.5 

20 

2.  187 

3  0 

0 

19.  1 

20 

2.  104 

0  2 

10 

25.8 

16 

2.070 

2  1 

7 

4.2 

6 

2.037 

1  0 

13 

7.5 

6 

2.000 

1  2 

8 

11.4  1 

20 

11 

4.2  1 

14 

1.957 

0  1 

14 

2.  1 

2 

1.872 

2  2 

0 

44.7 

45 

1.823 

2  1 

10 

6.4 

8 

1.799 

0  2 

13 

16.5 

14 

1.753 

1  2 

11 

2.5 

2 

1.725 

1  3 

4 

3.5 

3 

1.696 

3  1 

5 

6.9  1 

20 

1. 664 

2  1 

14 

15.  1  1 
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J.  EFFECT  OF  PRESSURE  ON  CRYSTAL  STRUCTURE  OF  CsMnCl0, 

RbMnCl3,  AND  RbFeCl3  3 

The  structures  of  ABX^  compounds,  where  A  is  a  large  univalent  cation  with  1.0  <  r A  <  1.9  A 

and  B  is  a  smaller  divalent  cation  with  0.5  <  r 0  <  1.0  A,  can  be  described  in  terms  of  the  hex- 

&  2+  4+ 

agonal  or  cubic  close  packing  of  AX^  layers.  Among  the  ABF^  and  A  B  compounds  there 
are  structures  stable  at  atmospheric  pressure  in  which  the  packing  is  entirely  cubic  (perovskite 
structure)  or  entirely  hexagonal  (2H,  in  polytype  formalism),  but  there  are  also  numerous  exam¬ 
ples  of  three  structures  in  which  the  packing  is  partially  hexagonal  and  partially  cubic.  With  in¬ 
creasing  pressure,  these  tend  to  transform  into  structures  in  which  there  is  an  increasing  frac¬ 
tion  of  cubic  packing.3^' 

According  to  the  literature  on  ABCL  compounds,  all  the  atmospheric  pressure  forms  defi- 

5  41 

nitely  identified  have  either  the  perovskite  or  2H  structure  except  for  CsCdCl^,  which  has  the 

6H  structure  in  which  two-thirds  of  the  layers  are  cubic-close-packed  (stacking  sequence  cchcch). 

42 

However,  examination  of  the  data  for  CsMnCl^  and  RbMnCl^  indicated  that  the  cell  dimensions 
are  not  consistent  with  either  entirely  cubic  or  entirely  hexagonal  close  packing.  We  have  deter¬ 
mined  the  atmospheric  pressure  structures  of  these  two  compounds  and  have  investigated  the 

effect  of  pressure  on  their  structures  and  on  the  structure  of  RbFeCL,  which  is  2H  at  atmos- 

,  .  43  3 

phene  pressure. 

Samples  of  CsMnCl^  were  prepared  by  the  reaction  of  stoichiometric  quantities  of  CsCl  and 
MnCl^  •  411^0  in  a  stream  of  titanium-scrubbed  N2*  The  reactants  were  slowly  heated  to  200 °C 
and  held  there  a  few  hours  to  remove  all  the  water  of  hydration,  after  which  the  temperature  was 
raised  to  500  °C  for  a  few  more  hours.  The  same  procedure  was  used  to  prepare  RbMnCl^  and 
RbFeCl^  from  the  appropriate  reactants.  The  two  Mn  compounds  are  pink  and  the  Fe  compound 
brown.  None  is  noticeably  hygroscopic. 

The  X-ray  diffraction  data  for  the  CsMnCL  phase  stable  at  atmospheric  pressure,  which 

3  42 

are  in  good  agreement  with  the  published  values,  show  that  this  phase  has  the  rhombohedral 
9R  structure,  in  which  one-third  of  the  layers  are  cubic-close-packed  (stacking  sequence 
chhchhchh).  Table  II- 5  compares  the  observed  intensities  and  those  calculated  by  assuming  ideal 
close-packed  positions  for  all  the  atoms  in  the  9R  structure  and  an  overall  temperature  factor  of 
1.5.  The  agreement  is  good,  which  indicates  that  a  detailed  refinement  would  reveal  the  devia¬ 
tions  from  the  ideal  atom  positions  to  be  small,  as  they  are  for  the  9R  forms  of  CsCoF^  (Ref.  37) 
and  BaRuO^  (Ref.  44). 

The  effect  of  pressure  on  the  CsMnCL  structure  was  determined  by  the  same  method  used 
36  3 

in  previous  investigations.  In  a  typical  run,  a  gold  encapsulated  specimen  was  subjected  to  the 
desired  pressure  in  a  belt  apparatus  and  then  heated  to  the  desired  temperature  by  means  of  an 
internal  graphite  heater.  After  half  an  hour  at  this  temperature,  the  sample  was  quenched  to 
room  temperature  before  the  pressure  was  released.  X-ray  diffraction  measurements  were  then 
made  at  atmospheric  pressure  on  the  quenched  sample. 

When  pressure  is  applied  to  the  9R  form  of  CsMnCl^,  it  is  transformed  first  to  the  6H  poly¬ 
type  and  then  to  the  perovskite  form.  Since  the  perovskite  structure  is  entirely  cubic,  no  further 

35  36 

transformation  can  occur.  The  same  progression  is  observed  for  CsCoF^  and  BaRuO^.  The 
results  of  the  high  pressure  experiments  are  summarized  by  the  pressure-temperature  phase 
diagram  shown  in  Fig.  11-20,  which  shows  the  stability  range  for  each  of  the  three  polytypes.  The 
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PRESSURE  (kbor*) 


Fig.  11-20.  Pressure-temperature  phase 
diagram  for  CsMnCI^. 


transformation  pressures  (~  lOkbars  for  9R  ^  6H,  ~  25  kbars  for  6H  ^  perovskite),  like  those 
for  all  other  high  pressure  transformations  of  CsBF^  compounds,  are  not  very  sensitive  to 
temperature. 

Table  II-6  compares  the  observed  X-ray  intensities  for  the  6H  polytype  of  CsMnCl^  with 
those  calculated  by  assuming  ideal  close-packed  positions  for  all  atoms  and  an  overall  tempera¬ 
ture  factor  of  1.0.  Some  of  the  discrepancies  between  observed  and  calculated  values  are  out¬ 
side  experimental  error,  indicating  that  the  positions  of  some  of  the  atoms  are  significantly 
different  from  the  ideal  ones.  Structure  refinements  for  other  6H  polytypes,  including  those  of 
RbNiF^  (Refs.  37  and  40)  and  CsCdCl^  (Ref.  41)  have  revealed  significant  distortions  from  the 
ideal  positions. 

Table  II- 7  gives  the  X-ray  results  for  the  perovskite  form  of  CsMnCl^.  In  this  case  there 
is  good  agreement  between  the  observed  intensities  and  those  calculated  by  assuming  ideal  atom 
positions  and  an  overall  temperature  factor  of  1.5. 

Both  high  pressure  forms  of  CsMnCl^  are  moderately  stable  at  atmospheric  pressure.  There 
was  no  evidence  of  retransformation  to  the  9R  structure  after  annealing  in  for  3  days  at  200 aC, 
but  complete  retransformation  occurred  after  only  a  few  hours  at  300°C.  These  retransformation 
temperatures  are  similar  to  those  found  for  CsBF^  high  pressure  polytypes,  which  retransform 
between  200°  and  500°C.  This  is  in  contrast  to  the  more  stable  high  pressure  ABO^  polytypes, 
some  of  which  retransform  slowly  even  at  1100°C. 

The  x-ray  data  for  the  form  of  RbMnCB  stable  at  atmospheric  pressure,  which  are  in  good 

42  ^ 

agreement  with  the  published  values,  show  that  this  phase  has  the  6ll  structure.  The  x-ray 
pattern  is  very  similar  to  the  one  given  in  Table  II- 6  for  the  high-pressure  6H  polytype  of 
CsMnCly  At  7 00°C,  the  6H  form  of  RbMnCl^  is  transformed  to  the  perovskite  structure  when 
the  pressure  exceeds  7  kbars.  Table  II- 7  shows  that  the  observed  x-ray  intensities  are  in  good 
agreement  with  those  calculated  by  assuming  the  ideal  perovskite  atom  positions  and  an  overall 
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TABLE  11-6 

X-RAY  DATA  FOR  HIGH  PRESSURE  CsMnCI3(6H) 

c 

i=  7. 268(5)  A,  c=  17.  85(5)  A 

hkl 

^calc 

^obs 

^obs 

hkl 

^calc 

^obs 

^obs 

1  0  2 

1.6 

4 

5. 12 

30  3 

0.  1 

2  1  5 

1.0 

26 

1.984 

0  0  4 

1.1 

2 

4.43 

20  7 

25.0 

1  0  3 

30.8 

33 

4.  30 

30  4 

1  1  8 

1-2J 

7. 9  i 

8 

1.898 

1  1  0 

1  0  4 

100.0 

72.6 

175 

3.48 

2  2  0 
20  8 

54.1  | 
17.01 

1 

75 

1.819 

1  1  2 

2.3 

— 

— 

2  1  7 

18.8 

15 

1.743 

2  0  1 

1  0  5 

14.4  | 
1.3  I 

i 

17 

3.10 

2  2  4 

0.5 

3  1  3 

5.0 

15 

1.679 

20  2 

42.  3  | 

I 

70 

2.961 

20  9 

10.4 

0  0  6 

27.2  1 

1 

3  1  4 

14.4| 

1 

15 

1.625 

1  1  4 

3.9 

- 

- 

2  1  8 

2.4] 

t 

20  3 

86.7 

90 

2.  778 

40  1 

1.61 

1 

3  1  5 

0.4 

t 

1 

6 

1.569 

2  0  4 

94.9 

90 

2. 566 

10  11 

5.7l 

1 

2  1  1 

0.3  1 

40  2 

5.1 

20  5 

12.5 

32 

2.360 

2  2  6 

19.5 

22 

1.552 

1  0  7 

24.0  1 

2  0  10 

8.2 

2  1  2 

1  1  6 

0.4  | 
1.2] 

i 

4 

2.297 

30  8 

3.8 

5 

1.527 

0  0  8 

1.4 

- 

- 

40  3 

2  1  9 

11.9| 

1.0] 

1 

15 

1.522 

2  1  3 

12.9 

12 

2.209 

40  4 

0  0  12 

15.3| 

5.4] 

1 

18 

1.487 

300 

22.7 

2  1  4 

31.1 

58 

2. 099 

1  0  8 

3.1 
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TABLE  11-7 

X-RAY  DATA  FOR  THE  HIGH  PRESSURE  PEROVSKITE  FORMS  OF  CsMnCI3  AND  RbMnCI3 

CsMnCI 

3  (a  =  5.  111(3)  A) 

RbMnCI^ 

,  (a  =  5.058(3)  A) 

hkl 

*calc 

*obs 

^obs 

^calc 

*obs 

^obs 

1  0  0 

6 

8 

5.  11 

1.4 

1 

5.05 

1  1  0 

100 

100 

3.62 

84 

85 

3.57 

1  1  1 

65 

61 

2.95 

66 

70 

2.915 

20  0 

81 

82 

2.55 

100 

100 

2.526 

2  1  0 

4 

4 

2.285 

0.4 

0 

- 

2  1  1 

44 

44 

2.087 

37 

40 

2.067 

2  2  0 

51 

47 

1.805 

62 

55 

1. 791 

300 

0.41 

1.705 

0.01 

2  2  1 

1-61 

2 

0.  ll 

0 

3  1  0 

18 

19 

1.616 

15 

16 

1.595 

3  1  1 

22 

22 

1.541 

21 

16 

1.520 

2  2  2 

16 

16 

1.475 

19 

16 

1.457 

3  2  0 

0.9 

0 

- 

0.0 

0 

- 

3  2  1 

20 

18 

1.366 

16 

14 

1.350 

4  0  0 

7 

8 

1.278 

8 

9 

1.262 

4  1  0 

3  2  2 

0.61 

0.6| 

0 

- 

0  0 

0  0 

0 

- 

4  1  1 

3  3  0 

3} 

9 

1.205 

21 

8 

1.  193 

3  3  1 

8 

8 

1.  172 

8 

7 

1.  159 

420 

18 

16 

1.  143 

21 

16 

4  2  1 

0.9 

1 

1.117 

0.0 

0 

- 

3  3  2 

4 

4 

1.090 

3 

2 

4  2  2 

13 

10 

1.044 

16 

10 

5  0  0 

0.  1  1 

0.0) 

4  3  0 

0.4) 

0 

0.0) 

0 

5  1  0 

4  3  1 

in 

9 

1.003 

?| 

5 

6) 

5j 

5  1  1 

3  3  3 

?i 

5 

0.983 

i) 

2 
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temperature  factor  of  1.5.  The  diffraction  pattern  is  seen  to  be  very  similar  to  the  one  given 
in  Table  II-7  for  the  perovskite  form  of  CsMnCl^.  At  atmospheric  pressure  the  perovskite  form 

of  RbMnCl-  retransforms  to  the  6H  structure  when  annealed  for  a  few  hours  at  400°C. 

J  43 

The  stable  form  of  RbFeCl^  at  atmospheric  pressure  has  the  2H  structure,  and  should 

have  high  pressure  polytypes.  The  results  of  high  pressure  studies  at  7  00°C  are  summarized 
in  Table  II-8.  No  single-phase  pattern  could  be  obtained  for  the  9R  structure  due  to  its  very  nar¬ 
row  stability  range.  The  transformation  pressures  are  not  sensitive  to  temperature  between 
400°  and  700°C. 

Preliminary  results  show  that  when  the  atmospheric  pressure  forms  of  the  materials  dis¬ 
cussed  above  are  prepared  in  HC1  and  not  they  are  very  hygroscopic.  There  is  no  change 
in  their  x-ray  patterns  and  in  the  AMnCl^  compounds,  there  is  no  change  in  any  of  the  transfor¬ 
mation  pressures.  For  RbFeCl^,  however,  all  the  transformation  pressures  are  lowered  by 
about  lOkbar  and  the  6H  —  3C  transformation  cannot  be  retained  by  quenching  to  room  tempera¬ 
ture.  We  think  now  that  use  of  gas  during  preparation  leads  to  nitride  ions  replacing  chloride 

ions.  The  extent  of  replacement  is  probably  very  small  in  the  Mn  compounds  but  larger  in 

2t 

RbFeCl^  with  its  easily  oxidized  Fe  ion. 


TABLE  11-8 

CRYSTALLOGRAPHIC  DATA  FOR  THREE  ABCI3  COMPOUNDS 

Stability  Range 

Cell  Dimensions 

Volume  (A^) 

Structure 

(kbars) 

(A) 

per  Formula  Unit 

CsMnCI3 

9R 

Atmospheric  to  10 

7.  288  X  27.  44  (Ref.  42) 

137.9 

6H 

10  to  25 

7.  268  X  1 7.  85 

136.  1 

Perovskite 

Above  25 

5.  115 

133.8 

RbMnCI3 

6H 

Atmospheric  to  ~7 

7. 1 64  X  17.798  (Ref.  42) 

131.8 

Perovskite 

Above  7 

5.058 

129.4 

RbFeCI3 

2H 

Atmospheric  to  25 

7. 060  X  6. 020  (Ref.  43) 

129.9 

9R 

25  to  30 

7.050  X  26.2 

125.3 

6H 

30  to  45 

7.04X  17.33 

124.0 

Perovskite 

Above  45 

4.965 

122.4 

For  all  three  ABCl^  compounds  investigated,  the  proportion  of  cubic  packing  in  the  high 
pressure  polytypes  increases  with  increasing  pressure.  The  crystallographic  data  for  all  the 
polytypes  of  these  compounds  are  summarized  in  Table  II -8 . 

J.  A.  Kafalas 
J.  M.  Longo 
D.  A.  Batson 
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K.  VOLUME  COMPRESSION  OF  TiOx 

Titanium  monoxide,  TiOx,  has  a  cubic  rock  salt  structure  which  is  stable  over  an  unusually 
wide  range  (0.75  <  x  <  1 . 30  at  1300°C)  and  has  an  unusually  high  vacancy  content  (up  to  16  per¬ 
cent  of  the  ideal  lattice  sites).  As  x  increases  from  0.75  to  1.30,  the  percentage  of  vacant  oxy¬ 
gen  sites  decreases  from  24  percent  to  0  and  the  percentage  of  vacant  metal  sites  increases 
from  8  to  24  percent.  Thus  the  total  number  of  vacancies  decreases  from  16  to  12  percent  with 

increasing  oxygen  content.  Even  for  the  stoichiometric  composition  TiO^  Q,  15  percent  of  both 

4  5 

metal  and  oxygen  sites  are  vacant.  Some  of  the  vacancies  can  be  removed  by  annealing  TiO 

46  47  x 

at  pressures  from  50  to  lOOkbars  and  temperatures  of  1100°  to  1600°C.  '  Because  of  the 

large  and  variable  number  of  vacancies,  some  of  which  can  be  removed  by  pressure  anneal, 

TiO^  might  be  expected  to  have  a  substantial  and  composition-dependent  volume  compression. 

We  are  studying  the  volume  compression  of  TiO  at  various  values  of  x  by  using  several 

48  X 

diamond  anvil  high  pressure  x-ray  units  to  measure  the  lattice  parameter  as  a  function  of 
pressure.  Two  of  these  units  are  gas  loaded,  a  third  has  screw  loaded  flat  springs  and  the 
fourth  has  direct  screw  loading.  The  internal  pressure  is  determined  by  observing  the  decrease 

in  lattice  parameter  of  NaCl  which  is  mixed  with  the  TiO  .  The  pressure  value  is  obtained  from 
49  x 

the  tabulated  results  of  pressure- volume  calculations  for  NaCl  made  using  equations  of  state 
and  the  elastic  constants  at  one  atmosphere.  Some  of  the  problems  associated  with  this  type  of 
experiment  are  described  in  a  recent  review/*8 

A  sample  of  -100  mesh  TiO^  is  mixed  with  dry  -200  mesh  NaCl  in  the  ratio  2TiOx:lNaCl, 
and  the  mix  is  prepressed  into  a  wafer  ~ 0.005  to  0.010  inch  thick.  Small  pieces  of  the  wafer 
are  glued  to  the  fixed  or  ''anvil11  diamond  of  the  squeezer.  An  atmospheric  pressure  reference 
pattern  is  made  on  a  portion  of  the  wafer  by  the  standard  Debye-Scherrer  technique,  and  the 
lattice  constants  of  both  compounds  are  determined  by  using  a  modified  Cohen  extrapolation.  An 
atmospheric  pressure  pattern  is  also  attempted  in  the  diamond-anvil  squeezer.  However,  it  is 
not  possible  to  know  when  the  diamonds  are  just  touching  the  wafer  so  that  the  first  film  may  be 
made  at  a  pressure  of  0.5  to  5  kbars.  A  series  of  patterns  is  then  made  at  increasing  loads. 


Fig.  11-21 .  Lattice  parameter  of  TiO^  q  as  a  function  of  pressure. 
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estimated  to  give  5  to  20  kbars  pressure  increments,  until  the  working  limit  of  the  unit  is  reached 

(100  to  150kbars)  or  until  there  is  evidence  of  a  failure  in  the  sample  or  diamonds.  If  a  series 

is  completed  without  a  failure,  a  final  pattern  is  made  with  all  the  load  removed  to  check  that 

the  fixed  diamond  has  not  moved  in  relation  to  the  film  during  the  pressure  cycle. 

The  x-ray  films  for  samples  under  pressure  require  exposures  of  100  to  400  hours  to  Zr- 

filtered  Mo  radiation  from  a  Norelco  fine  focus  x-ray  tube,  or  exposures  of  15  to  40  hours  to 

radiation  from  a  Jarrell-Ash  microfocus  unit  (used  through  the  courtesy  of  the  Jarrell-Ash  Co.). 

The  better  films  give  four  to  six  diffraction  lines  for  each  compound  and  the  best  films,  which 

50 

are  obtained  in  the  Bassett  screw  loaded  squeezer,  have  seven  or  eight  lines  for  each  with 
29  as  high  as  44°.  For  each  compound  a  separate  lattice  parameter  value  is  calculated  from 
the  d-value  of  each  line.  The  lattice  constant  adopted  is  usually  the  numerical  average  of  all 
the  individual  values.  However,  on  occasion  one  or  more  low  angle  lines  may  appear  grossly 
inconsistent  and  are  not  used. 

Data  for  a  typical  series  of  measurements  are  shown  in  Fig.  11-21,  where  the  lattice  con¬ 
stant  of  TiO^  q  is  plotted  vs  pressure  up  to  108  kbars.  This  series  was  made  in  two  different 
gas  loaded  squeezers  using  portions  of  the  same  TiO^  ^-NaCl  mix.  The  values  of  volume  com¬ 
pression  (AV/VoAp)  for  this  and  the  other  compositions  studied  are  listed  in  Table  II-9.  Each 
value  is  based  on  the  best  straight  line  obtained  by  a  visual  fit  to  all  the  data  points  for  that 
composition  as  shown  in  Fig.  11-21,  even  though  in  some  cases  an  individual  run  could  be  better 

described  by  a  curve  indicating  a  decrease  in  compression  at  higher  pressures.  The  differences 

-4  -4  -1 

among  the  tabulated  values,  which  range  from  2.6  X  10  to  4.0  x  10  kbar  ,  are  probably  not 
significant  in  view  of  the  experimental  errors. 


TABLE  11-9 

VOLUME  COMPRESSION  OF 

TiO 

X 

X 

Pressure 

Range 

(kbars) 

No.  of  Runs 

-AV 

V6Ap 

(kbar'1) 

0.83 

0-124 

2 

3. 3  X  10'4 

0.85 

0-39 

3 

2. 3  X  lO-4 

1.00 

0-108 

2 

3. 8  X  10'4 

1.10 

0-139 

1 

2.6  X  10"4 

1.25 

0-52.5 

1 

4.  OX  10'4 

There  are  two  other  reports  of  volume  compression  measurements  on  TiOx,  both  carried 

out  with  tetrahedral-anvil  pressure  units,  which  generate  nearly  hydrostatic  pressure.  Taylor 

and  Doyle51  report  a  value  of  5.5(8)  x  10“4kbar_1  for  TiO.  n  (with  15  percent  vacancies)  with- 

52  1 ' U 

out  giving  any  experimental  details.  Iwasaki,  et  al.,  measured  the  lattice  parameter  for 
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TiOi  25  to  60kbars  using  NH4CI  as  a  pressure  marker.  They  give  a  compression  value  of 

8.6  x  lO^kbar"1  based  on  their  data  above  lOkbars.  (Between  1  atmosphere  and  1  0  kbars  they 

observed  a  slight  increase  in  lattice  parameter,  but  this  is  such  an  improbable  result  that  all 

possible  sources  of  experimental  error  must  be  eliminated  before  it  can  be  accepted.  Iwasaki, 
52 

et  al. ,  do  not  report  the  diffraction  lines  they  observed  or  give  any  analysis  of  errors.) 

Although  the  quantitative  values  from  the  three  investigations  of  TiOx  are  somewhat  differ¬ 
ent,  it  is  clear  that  in  spite  of  the  large  number  of  vacancies  the  volume  compression  of  TiOx 
is  very  small  and  not  markedly  dependent  on  the  type  or  number  of  vacancies.  In  fact,  its  com¬ 
pression  is  close  to  those  for  TiC  and  TiN  and  of  the  same  order  as  those  for  monoxides  such 
as  MnO  and  CoO  which  have  very  few  vacancies.  M  D  Banus 

Mary  C.  Lavine 


REFERENCES 

1.  T.B.  Reed  and  E.R.  Pollard,  ].  Cryst.  Growth  2,  243  (1968),  DDC  AD-681586. 

2.  G.M.  Grover,  T.P.  Cotter  and  G.  F.  Erikson,  ].  Appl.  Phys.  35,  1990  (1964). 

3.  G.Y.  Eastman,  Sci.  Amer.  218,  28  (1968). 

4.  F.P.  Koffyberg,  Phys.  Letters  A30,  37  (1969)  and  private  communication. 

5.  M.  Altwein,  H.  Finkenrath,  C.  Konak,  ].  Stuke  and  G.  Zimmerer,  Phys.  Stat.  Solidi  29, 

203  (1968). 

6.  S.  Hayashi,  Oyo  Butsuri  34,  825  (1965). 

7.  P.  Hoschl,  C.  Konak,  and  V.  Prosser,  Mat.  Res.  Bull.  £,  87  (1969). 

8.  Solid  State  Research  Report,  Lincoln  Laboratory,  M.I.T.  (1969:3),  p.  15,  DDC  AD-696620. 

9.  Solid  State  Research  Report,  Lincoln  Laboratory,  M.I.T.  (1969:4),  p.  8,  DDC  AD-701022. 

10.  R.  Hultgren,  R.  L.  Orr,  P. D.  Anderson  and  K.  K.  Kelley,  Supplement  to  Selected  Values  of 
Thermodynamic  Properties  of  Metals  and  Alloys  (Wiley,  New  York,  1963). 

11.  J.M.  Haschke  and  H.A.  Eick,  J.  Phys.  Chem.  73,  374  (1969). 

12.  C.E.  Holley,  E.J.  Huber  and  F.  B.  Baker  in  Progress  in  the  Science  and  Technology  of  the 
Rare  Earths  3,  edited  by  L.  Eyring  (Pergamon  Press,  Oxford,  1968),  p.  343. 

13.  S.  Larach,  R.E.  Shrader  and  C.  F.  Stocker,  Phys.  Rev.  108,  587  (1957). 

14.  R.E.  Halsted,  D.T.F.  Marple  and  M.  Aven,  Extended  Abstracts,  Electronics  Div., 
Electrochem.  Soc.  F4,  6  (1965). 

15.  A.C.  Aten,  J.  Phys.  Chem.  Solids  28,  1340  (1967). 

16.  G.W.  Iseler  and  A.J.  Strauss,  J.  Luminescence  3,  1  (1970). 

17.  M.  Cardona,  K.  L.  Shaklee  and  F.H.  Poliak,  Phys.  Rev.  154,  696  (1967). 

18.  A.G.  Thompson  and  J.C.  Woolley,  Canad.  J.  Phys.  45,  255  (1967). 

19.  J.A.  Van  Vechten  and  T.  K.  Bergstresser ,  Bull.  Am.  Phys.  Soc.  15,  254  (1970). 

20.  L.R.  Shiozawa  and  J.M.  Jost,  "Research  on  II  —VI  Compound  Semiconductors,"  Clevite 
Corporation  Electronic  Research  Division,  Cleveland,  Ohio  (May  1965),  p.  118, 

DDC  AD -620297. 


44 


Section  II 


21.  Solid  State  Research  Report,  Lincoln  Laboratory,  M.I.T.  (1969:3),  p.  21,  DDC  AD-696620. 

22.  J.  Steininger,  Mat.  Res.  Bull.  3,  595  (1968),  DDC  AD-674790. 

23.  D.T.F.  Marple,  Phys.  Rev.  150,  728  (1966). 

24.  O.  Brafman  and  I.T.  Steinberger,  Phys.  Rev.  143,  501  (1966). 

25.  L.  V.  Prytkina,  V.  V.V.  Voiko,  A.J.  Mentser,  A.V.  Vanyukov  and  P.  S.  Kireev,  Sov.  Phys.  - 
Semicon.  2,  509  (1968). 

26.  D.B.  Rogers,  R.J.  Arnott,  A.  Wold  and  J.  B.  Goodenough,  J.  Phys.  Chem.  Solids  24,  347 
(1963). 

27.  D.B.  Rogers,  J.  L.  Gillson  and  J.E.  Gier,  Solid  State  Comm.  5,  263  (1967). 

28.  B.  Reuter  and  K.  Muller,  Naturwiss.  7,  164  (1967). 

29.  J.  Jaskowsky,  Ber.  Bunsengesell.  Phys.  Chem.  70,  199  (1966). 

30.  J.  B.  Goodenough,  J.  Solid  State  Chem.  (in  press). 

31.  W.  Riidorff,  J.  Kandler  and  D.  Babel,  Z.  anorg.  allgem.  Chem.  317,  261  (1962). 

32.  R.D.  Burbank  and  H.T.  Evans,  Acta  Cryst.  1_,  330  (1948). 

33.  J.  A.  Kafalas  and  J.  M.  Longo,  Mat.  Res.  Bull.  3,  501  (1968),  DDC  AD-674785. 

34.  Y.  Syono,  S.  Akimoto  and  K.  Kohn,  J.  Phys.  Soc.  Japan  26,  993  (1969). 

35.  J.  M.  Longo  and  J.  A.  Kafalas,  J.  Solid  State  Chem.  1,  103  (1969),  DDC  AD-693900. 

36.  _ ,  Mat.  Res.  Bull.  3,  687  (1968),  DDC  AD-674791. 

37.  D.  Babel,  Z.  anorg.  allgem.  Chem.  369,  117  (1969). 

38.  Z.Zalkin,  K.  Lee  and  D.  Templeton,  J.  Chem.  Phys.  37,  697  (1962). 

39.  G.A.  Smolenskii,  MJ\  Petrov,  V.V.  Moskalev,  V. S.  L'vov,  V.S.  Kasperovich  and 
E.V.  Zhirnova,  Sov.  Phys.  Solid  State  10,  1040  (1968). 

40.  R.J.  Arnott  and  J.M.  Longo,  J.  Solid  State  Chem.  (in  press). 

41.  S.  Siegel  and  E.  Gebert,  Acta  Cryst.  J_7,  790  (1964). 

42.  M.  Kestigian,  W.J.  Croft  and  F.D.  Leipziger,  J.  Chem.  Eng.  Data  12,  97  (1967). 

43.  H.J.  Seifert  and  K.  Klatyk,  Z.  anorg.  allgem.  Chem.  342,  1  (1966). 

44.  P.  C.  Donohue,  L.  Katz  and  R.  Ward,  Inorg.  Chem.  4,  306  (1965). 

45.  M.D.  Banus  and  T.  B.  Reed,  Chemistry  of  Extended  Defects  in  Non-Metallic  Solids, 
edited  by  L.  Eyring  and  M.  O'Keefe.  (North-Holland,  Amsterdam,  1970),  p.  488. 

46.  M.D.  Banus,  Mat.  Res.  Bull.  3,  723  (1968),  DDC  AD-677839. 

47.  N.J.  Doyle,  J.K.  Hulm,  C.K.  Jones,  R.C.  Miller  and  A.  Taylor,  Phys.  Letters  26A,  604 
(1968). 

48.  M.D.  Banus,  High  Temperatures  —  High  Pressures  1^,  No.  5  (in  press). 

49.  D.L.  Decker,  W.A.  Bassett,  L.  Merrill,  H.T.  Hall  and  J.D.  Barnett,  "High  Pressure 
Calibration  -  A  Critical  Review,"  High  Pressure  Data  Center,  Brigham  Young  University, 
(October,  1969). 

50.  W.A.  Bassett,  T.Takahashi  and  P.W.  Stook,  Rev.  Sci.  Instr.  38,  37  (1967). 

51.  A.  Taylor  and  N.J.  Doyle,  Chemistry  of  Extended  Defects  in  Non-Metallic  Solids, 
edited  by  L.  Eyring  and  M.  O'Keeffe  (North-Holland,  Amsterdam,  1970),  p.  523. 

52.  H.  Iwasaki,  K.  Kamigaki,  S.  Ogawa,  O.  Terasaki  and  D.  Watanabe,  J.  Phys.  Soc.  Japan 
(in  press). 


45 


III.  PHYSICS  OF  SOLIDS 


A.  ELECTRONIC  BAND  STRUCTURE 

1.  Fermi  Surface  and  Optical  Properties  of  Potassium 

We  have  carried  out  a  new  calculation  of  the  interband  absorption  of  potassium.  Until  now 

1 

the  optical  constants  have  been  calculated  using  a  local  pseudopotential  model.  However,  it 
2 

has  been  demonstrated  that  the  local  pseudopotential  formulation  is  inadequate  to  account  for 
the  experimentally  measured  Fermi  surface  of  potassium.  It  has  been  shown  that  the  pseudo¬ 
potential  that  appears  in  the  optical  matrix  element  will  incorporate  the  effects  of  many-body 
correlations  only  if  it  is  consistent  with  the  experimental  data  of  Fermi  surface  measurements 
(such  as  dHvA  results).  The  two  nonlocal  pseudopotentials  used  here  were  selected  with  this 
requirement  in  mind. 

One  of  the  nonlocal  pseudopotentials  V-*(G  —  G’)  used  (where  G  is  a  reciprocal  lattice  vector 
—  k 

and  k  the  electron  momentum)  is  a  simple  generalization  of  the  local  pseudopotential,  consistent 
with  bcc  symmetry.  It  is  parameterized  phenomenologically  by  constants  characterizing  its 
local  and  nonlocal  (k  ■  G-dependent)  features.  The  other  potential  <k  +  G|r|k  +  Gf>  is  the 

4 

Ziman  form  of  the  KKR  method  in  which  the  wavefunction  external  to  the  muffin- tin  radius  R 

s 

is  expressed  as  plane- waves.  It  is  parameterized  by  s-,  p- ,  and  d- phase  shifts.  R  was  taken 
to  be  the  atomic  core  radius. 

Each  potential,  when  inserted  into  a  19  X  19  determinant  and  solved  for  the  energy  eigenval¬ 
ues,  gave  an  excellent  fit  (Fig.  Ill- 1 )  to  the  experimental  Fermi  surface.  The  KKR-Z  result 
indicated  that  the  d-phase  shift  was  particularly  important  —  which  is  probably  the  reason  for 

the  nonlocality  required  of  any  pseudopotential  used  in  the  case  of  potassium.  It  is  of  interest 

5 

that  each  potential  was  weaker  than  the  Ashcroft  local  pseudopotential  and  the  electronic  states 
are  thus  more  free-electron-like. 

While  permitting  a  good  fit  to  the  Fermi  surface,  nonlocal  pseudopotentials  suffer  a  disad¬ 
vantage  over  local  potentials  when  it  comes  to  calculating  properties  dependent  on  the  wavefunc- 
tions.  In  general,  the  wavefunctions  are  no  longer  expressed  in  terms  of  plane  waves.  However, 
the  Ziman  form  of  the  KKR  method  is  constructed  so  that  the  eigenfunctions  are  expressible 

simply  as  plane  waves  outside  R  .  In  the  case  of  V_*(G  —  G’),  this  is  not  necessarily  true. 

s  k 

The  optical  conductivity  <j(o;)  was  calculated  by  solving  the  19  x  19  determinant  at  different 
points  in  the  Brillouin  zone  as  required  by  the  Monte  Carlo  integration  scheme  used.  With  either 
nonlocal  pseudopotential,  the  interband  threshold  occurs  at  about  1.3  eV  and  the  absorption  max¬ 
imum  at  about  2  eV  (Fig.  III-2).  The  absorption  magnitude  depends  on  the  potential,  however 
13  - 1 

Re  o( 2  eV)  ~  1.8  X  10  sec  with  the  potential  V,  and  five  times  smaller  with  the  KKR-Z  poten¬ 
tial.  The  former  result  is  closer  to  the  measurements  of  Smith.^  The  source  of  the  discrepancy 
seems  most  likely  to  occur  in  the  KKR-Z  calculation  on  account  of  the  approximations  made  in 
the  calculation  —  namely,  neglect  of  the  energy  dependence  of  the  pseudopotentials  and  truncation 
of  the  eigenfunctions  to  19  plane  waves.  Consideration  of  these  features  should  be  included  in 
further  calculations.  A.  R.  WilSon 

C.  Y.  Young 
G.  F.  Dress  el  haus 
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ANGLE  (d«g) 

Fig.  Ill-l.  Best  fit  to  experimental  fractional  change  AA/A  of  central  cross  sectional 
area  with  direction. 


Fig.  MI-2.  Interband  absorption,  Re  ct(u),  for  the  two  pseudopotentials  used  and  the  experimental 
point  of  N.V.  Smith  (Ref.6). 
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2.  Model  of  the  Insulator-Metal  Transition  in  TigOg 

In  an  attempt  to  understand  the  nature  of  the  insulator- metal  transition  in  Ti^O^  (Ref.  7),  we 
have  explored  a  model  based  on  the  competition  between  elastic  and  electronic  free  energies  of  the 
system  as  a  function  of  temperature.  The  model  consists  of  two  bands,  one  filled  and  one  empty 
at  T  =  0°K,  separated  by  an  energy  gap  which  depends  linearly  on  a  lattice  distortion  parameter. 
The  normalized  free  energy  of  the  system  is  written. 


F  =  62  —  *jkT  ^  In  [1  +  e"^€“^)AT]  PeU)  dc 


+  kT  ^  In  [1  +  e"(e+^1"6l,/kT]  ph(e)  dc  j 


(III- 1 ) 


In  Eq.  ( III—  1 ),  6  is  a  normalized  variable,  proportional  to  the  change  in  lattice  distortion  param¬ 
eter  at  T  =  0°K.  The  first  term  is  the  lattice  distortion  free  energy,  the  second  is  the  free 
energy  due  to  the  excitation  of  electrons  into  the  empty  band, 
and  the  third  term  is  the  free  energy  of  the  holes  excited  in  the 
filled  band  (see  Fig.  Ill- 3).  e  is  a  normalized  energy  variable, 

(jl,  the  normalized  Fermi  energy,  and  [1  —  6],  the  normalized 
energy  gap.  All  energies  and  kT  have  been  normalized  to  the 
energy  gap  at  T  =  0°K.  Pe<0  and  P^U)  are  the  densities  of 
states  for  the  bands  which  are  empty  and  filled,  respectively, 
at  T  =  0°K.  is  a  dimensionless  parameter  which  determines 
the  relative  scale  of  the  lattice  distortion  and  electronic  free 
energies. 

With  this  expression  for  the  free  energy,  both  6  and  \i  can 
be  determined  at  any  temperature  by  the  conditions  that  8F/86  = 

0;  and  the  number  of  electrons  equals  the  number  of  holes. 

Such  a  model  is  capable  of  giving  a  reasonable  description 
of  the  change  in  c/a  ratio  in  Ti^O^  as  a  function  of  temperature. 

As  an  example,  we  show  in  Fig.  Ill- 4  the  results  of  a  computer 
calculation  of  6  vs  kT. 

The  same  model  should  also  describe  the  behavior  of  a 
Raman  mode  frequency  associated  with  the  displacement  6, 
as  a  function  of  temperature,  through  the  relation,  CJj^aman  00 


Fig.  MI-3.  Schematic  population 
of  band  model  at  T  =  0°K. 


Jd2F/i 


d6 


2  2 

A  plot  of  d  F/d6  ,  shown  in  Fig.  Ill- 5,  predicts  a  very  pronounced  "mode  softening" 


in  the  vicinity  of  the  transition,  due  to  the  strong  coupling  of  the  electronic  and  lattice  distortion 
free  energy  needed  to  explain  the  sharp  change  in  c/a  ratio  observed  in  T^O^  in  the  transition 
region  since  the  Raman  frequencies  observed  in  Ti^O^  (Ref.  8)  show  only  a  slight  softening,  we 
can  conclude  that  this  model  based  on  lattice  distortion  as  the  driving  mechanism  for  the  transi¬ 
tion  in  Ti^O^  is  inconsistent  with  the  Raman  scattering  data. 


P.  M.  Raccah 
H.  J.  Zeiger 
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Fig.  111-4.  Theory  curve  shows  8  vs  kT  for  normalized  upper  bandwidth  =  5. 0,  lower 
bandwidth  =  0.2  and  k  =  37.  Data  curve  shows  c/a  behavior  as  function  of  temperature 
on  same  normalized  scale. 


Fig.  II 1-5.  Comparison  between  experimental  phonon  frequencies  and  theoretical  predictions. 
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3.  Variational  Approach  to  the  Metal -Semiconductor  Transition 

A  mechanism  for  a  metal  to  semiconductor  transition  occurring  with  increasing  temperature 

9 

has  been  found  in  the  context  of  a  new  variational  approximation  to  the  half- filled  band  Hubbard 

Hamiltonian.  We  have  found  the  stability  boundary  (based  on  minimum  free  energy)  between  the 

intuitively  simple  thermodynamic  states  representing  localized  (Wannier  state)  and  band  (Bloch 

state)  behavior.  This  is  shown  in  Fig.  Ill- 6,  in  which  I  is  the  Coulomb  repulsion,  A  is  a  meas- 

2  2 

ure  of  the  bandwidth  defined  by  A  =  2  b. .  is  the  transfer  integral  between  sites  i  and  j  and 


Ha-8-mHl 


Fig.  III-6.  Stability  boundary  between 
Wannier  (localized)  and  Block  (band) 
states. 


WANNIER 


BLOCH 


1 0 

bii  =  0.  While  this  boundary  had  previously  been  given  for  a  linear  chain  with  nearest- 
neighbor  b. ^  the  present  boundary  was  calculated  using  the  density  of  states 


P<€)  *  2  A  ^  UI<2A 

=  0  |  e  |  >  2A  (III- 2) 

1 1 

which  was  used  by  Hubbard  in  his  investigation  of  an  insulator- metal  transition  at  T  =  0°K. 
Equation  ( III- 2 )  is  not  unreasonable  for  a  more  realistic  3-dimensional  band  structure. 

We  make  the  following  points. 

1 0 

(a)  The  qualitative  behavior  is  very  similar  to  that  for  the  linear  chain, 
and  indeed  can  be  understood  in  general  terms  as  follows: 

The  fact  that  the  low  temperature  slope  is  positive  can  be  understood 
easily  if  we  note  that  (1)  the  zero- temperature  entropy  is  Nk  In  2  for  the 
Wannier  state  and  zero  for  the  Bloch  state  and  (2)  the  ground- state 
energy  difference  between  the  Bloch  and  Wannier  states,  Ep  —  E^, 
decreases  with  increasing  bandwidth.  Thus  at  low  T,  the  free  energies 

=  E^y  —  kT  In  2  +.  .  . 


Fb  =  Eb  +  °(T2) 


(III-  3) 


so  that  on  the  boundary,  kT  In 2  =  E^-  E^.  Hence  on  the  boundary 


8T 

0A 


1 

k  In  2 


9(EB~  EW> 

8A 


>  0 


(HI- 4) 
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Also,  the  fact  that  the  boundary  asymptotes  to  a  vertical  line  can  be 
understood  when  it  is  considered  that  (1)  as  T  -*■  00  the  entropy  for  the 
Wannier  and  the  Bloch  states  are  the  same  (Nk  In 4)  and  (2)  the  inter¬ 
nal  energies  approach  a  constant  independent  of  T.  The  equality  of 
the  entropies  comes  from  the  fact  that  for  kT  »  U,  polar  states  are 
highly  probable  in  the  Wannier  state  as  well  as  in  the  Bloch  state. 

Clearly  the  positive  slope  implies  a  transition  from  band  to  localized 
behavior  with  increasing  T  for 


1 

4a 


(b) 


(c) 


where  (A/I)q  is  the  zero- temperature  intercept  of  the  boundary. 

The  phase  boundary  obtained  here  is  in  qualitative  disagreement  with 
one  presented  by  Mott 12  (the  slope  of  his  curve  is  negative).  The 
reason  for  the  difference  is  that  Mott  considered  a  transition  from  a 
magnetically- ordered  insulating  state  to  a  paramagnetic  metallic 
state  whereas  we  have  considered  paramagnetic  states  on  both  sides 
of  the  transition.  It  seems  to  us  that  in  a  consistent  treatment,  one 
should  consider  magnetically  ordered  states  on  both  sides  of  the 
transition  at  low  temperatures,  as  we  have  done.10  The  present 
treatment  is  restricted  to  temperatures  greater  than  any  magnetic- 
ordering  temperature. 

We  note  that  the  zero-temperature  intercept  (A/l)0  =  0.295,  is  very 
close  to  the  value  0.289  at  which  Hubbard11  obtained  a  nonmagnetic 
insulator-metal  transition  at  0°K. 

T.  A.  Kaplan 

R.  A.  Bari 


(HI- 5) 


4.  Interaction  of  Bound  Electrons  with  Local  and  Resonant  Modes 
in  Semiconductors 

It  is  suggested  that  some  of  the  multiple  pinning  phenomena  observed  in  the  optical  and 

13  14 

infrared  spectra  of  impurities  and  other  imperfections  in  semiconductors  ’  may  be  due  to 
pinning  to  local  or  resonant  modes.  As  a  model  problem,  we  have  considered  the  optical  ab¬ 
sorption  spectrum  of  an  electron  in  an  insulator  bound  to  an  impurity,  interacting  with  the  L-O 
branch  of  the  lattice  vibrational  spectrum,  as  affected  by  the  presence  of  the  impurity.  A 
Green's  function  calculation  for  the  spectrum  in  the  case  of  a  local  mode,  gives  a  pinning  effect 
due  to  interaction  with  the  local  mode  in  the  vicinity  of  the  impurity,  as  well  as  pinning  due  to 
the  lattice  modes.  The  nature  of  the  pinning  phenomena  in  the  case  of  a  resonant  mode,  and  the 
criteria  for  the  observation  of  such  pinning,  are  being  investigated. 

R.  W.  Davies 
H.  J.  Zeiger 


5.  Band  Theory  of  Optical  Activity 

Optical  activity  arises  from  terms  linear  in  q  in  the  wavevector-  and  frequency- dependent 
dielectric  tensor,  e(q,  co ).  In  crystals  c(q,  a> )  is  given  compactly  to  all  powers  of  q  (all  multipole 
moments)  using  a  Bloch  state  band  formalism.  The  linear  q  part  of  e  is  obtained  by  a  first  order 
q  •  p  perturbation  expansion  of  either  the  band  state's  energy  or  wavefunction.  For  the  uniaxial 
crystal  classes,  D^,  D^,  D^,  Cy  C^,  C^,  the  dielectric  tensor  component  responsible  for  the 
optical  activity  for  propagation  parallel  to  the  c-axis  is 
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c  (q  ,  cj  ) 
xy  ^z 


where 

47re2f,  (ky  |p  | k-y ' )  (ky1  |p  Iky) 

«i  s  “dr  - Vt - ""•6b’ 

m  fifi  -  a;  ) 

and  for  propagation  in  the  basal  plane 


y  4qz(ky' |pz|ky') 


kyy ' 
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y  y 


j  2  2 

3co  ,  —  co 

y  y 
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CO  ,  —  CO 

yy 


(III- 6a) 


e  (q  a) )  =  —  i  y.  2q  a[cM  —  cMl 
yzVMz  ^  Mx  1  e  o J 


(III-  7a ) 


where 

a  =  (ky'  | x  | ky " )  =  i(ky'  |px  |ky")/mu^,^,, 


{ III—  7b) 


4neZ{ky  (kyIPy  Ikv')  (ky"  |pz|ky) 


m2fi  n 


/  "2 - 27 

<ylUy«y-  W  ) 


(III—  7c ) 


and  is  given  by  an  expression  similar  to  Eq.  (Ill- 7c)  with  vyUy  w  i  •  Here  ^y\y  is  the 

energy  difference  between  Bloch  states  |ky’)  and  |ky);  and  f,  is  the  Fermi  distribution  function, 
- 1 

[exp(E^—  E^J/kT  +  1]”  .  The  rotary  powers  parallel  and  perpendicular  to  the  c-axis  are 


p  i ,  =  -  iq  € 
r  1 1  Mz  xy'  o 

p,  =  — iq  c  /2(c  e  )* ^ 
^x  yz'  o  e' 


(III- 8a) 
(III-  8b ) 


where  cq  and  are  the  ordinary  and  extraordinary  q  =  0  dielectric  constants.  These  expressions 
may  be  used  to  calculate  the  optical  activity  from  a  given  band  model  or,  conversely,  to  adjust 
the  band  parameters  to  fit  the  rotary  power  data. 

The  general  results  above  take  on  a  revealing  form  in  a  simple  four- band  model  of  sym¬ 
metry  with  a  direct  gap  at  the  T- point.  Here  the  filled  band  is  assumed  to  have  T ^  symmetry 
and  the  conduction  band,  to  consist  of  a  T 2  singlet  (transforming  as  z)  and  a  doublet  (trans¬ 
forming  as  x,  y).  Then  =  cq;  and,  if  we  also  assume  that  (T^  |p^  1 T^)  -  (T^  |pz  |r^), 

Zk  —  €"]  =  —  cq.  This  latter  "cubic"  approximation  is  equivalent  to  attributing  the  total 

birefringence  to  the  energy  denominators  and  not  to  the  matrix  elements.  The  spatial  matrix 
element,  a,  has  an  upper  limit  of  the  size  of  the  unit  cell.  Thus  the  pairwise  near- cancellation 
of  two  resonances,  manifest  in  the  birefringence,  yields  the  correct  order-of- magnitude  of  the 
optical  activity.  In  the  case  of  small  birefringence,  i.e.,  |cuy  |  «  ~  uy"y>  anisotropy 

of  the  rotary  power  is 


ill 

pi 


-2 


<r,|pzlr3) 

(r2lPx'r3> 


( III- 9 ) 
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TEMPERATURE  (K) 


Fig.  1 1 1  —7 (a ) .  Temperature  of  kink  in  magnetization  vs  temperature  curves  of  EuO 
as  a  function  of  applied  field  H  . 


tc-t  CK) 


Fig.  Ill-7(b).  Deviation  of  kink  temperature  from  T  as  a  function  of  applied  field. 
Slope  corresponds  to  (3  =  0.368. 
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If  the  "cubic"  approximation  is  again  invoked,  the  matrix  element  ratio  is  ~  1  since  + 
r3  =  ■  T  (triplet).  Then  p^  /p^  ~  —  2,  which  is  close  to  the  measured  ratio  in  a-quartz. 


A.  S.  Pine 

G.  F.  Dresselhaus 


B.  MAGNETISM 

1.  Critical  Magnetic  Properties  of  EuO 

A  study  is  being  made  of  the  magnetic  properties  of  single  crystal  EuOt  in  the  vicinity  of 
the  Curie  temperature.  We  have  analyzed  the  measurements  to  determine  the  magnetic  critical 
indices  p,  6  and  y,  which  are  defined  by  the  equations 


(7  <=  (T  -  T  f 
c 

for  T  <  T 

c 

(III- 9a) 

H  o6 

for  T  =  T 

c 

(III- 9b) 

X~ 1  «  (T—  TCP 

for  T  >  T 

c 

(III- 9c) 

where  <r  represents  the  magnetic  moment  of  the  sample,  H  is  the  magnetic  field  internal  to  the 
sample  (applied  field  minus  demagnetization  effects),  \Q  is  the  initial  susceptibility  and  Tc  is 
the  Curie  temperature. 

1 5 

The  critical  index  p  was  determined  by  the  method  proposed  by  Rayl  and  Wojtowicz, 
which  utilizes  the  fact  that  there  is  a  sharp  kink  in  the  magnetization  vs  temperature  curve  at 
the  temperature  where  the  applied  field  (H&)  is  equal  to  the  demagnetizing  field  of  the  sample. 

At  this  point  the  applied  field  is  proportional  to  the  magnetization  according  to  the  relationship 
Ha  =  NM,  where  N  is  the  demagnetizing  factor  (4tt/3  for  a  sphere).  This  kink  was  extremely 
sharp  in  our  crystal,  occurring  over  less  than  0.01  °K,  indicating  high  sample  homogeneity. 

The  locus  of  kink  points  as  a  function  of  applied  field  is  shown  in  Fig.  III-7a.  It  is  seen  that  ex¬ 
trapolation  of  this  curve  to  H  =0  gives  an  accurate  measure  of  T  .  Replotting  the  points  in 

a  c 

Fig.  III-7a  on  a  log- log  scale  gives  the  linear  relationship  predicted  by  Eq.  (Ill- 9a)  with  p  =  0.368, 

as  shown  in  Fig.  Ill- 7b .  The  relationship  is  seen  to  be  valid  over  at  least  two  decades,  from 

T  -  T  =  0.1  0  to  10°. 
c 

The  critical  temperature  exponent  6  is  determined  directly  from  measurement  of  the  mag¬ 
netic  moment  vs  internal  field  at  T  =  T  .  The  results  of  this  measurement,  as  given  in  Fig.  Ill- 8, 
correspond  to  6  =  4.46. 

The  initial  susceptibility  in  the  region  immediately  above  the  Curie  temperature  was  deter¬ 
mined  by  measuring  the  magnetic  moment  as  a  function  of  temperature  at  several  fields  and  ex- 
2 

trapolating  the  cr  vs  H/a  isotherms  obtained  to  H  =  0.  The  resultant  values  of  xn  were  analyzed 

o  |  ^ 

by  two  distinct  methods.  The  first  was  the  T*  technique  introduced  by  Kouvel  and  Fisher. 

This  involves  the  determination  of  T*  =  x”  */dx~  */dT)  as  a  function  of  temperature.  Defining  a 
temperature-dependent  effective  exponent  y*  (T)  =  (T  —  T  )/T*  ,  then  y*  (T)  =  y  in  the  limit 
T  =  T  .  A  computer  program  was  written  to  evaluate  y*  (T)  according  to  this  procedure,  and 


tThis  crystal  was  grown  by  T.  B.  Reed. 
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Fig.  111-8.  Magnetic  moment  a  vs  internal  magnetic  field  in  EuOat  T  =  T  . 
Slope  corresponds  to  6  =  4.46.  c 


Fig.  111-9.  Inverse  of  initial  susceptibility  x  vs  T  -  T  ,  assuming  T  =69.594°K. 
Slope  corresponds  to  y  =  1 .285. 


56 


Section  til 


it  was  found  that  for  Tc  =  69.59  °K,  the  value  of  y  was  essentially  constant  and  equal  to 
1.278  ±  0.002  for  T  —  Tc  between  0.6  and  4°K.  At  higher  temperatures  y*  decreases  with  in¬ 
creasing  temperature. 

.  A  second  method  involved  the  assumption  that  y  remains  constant  over  the  temperature 
interval  being  measured,  and  the  use  of  a  computer  to  obtain  a  direct  determination  of  the  best 
values  of  y  and  Tc  on  the  basis  of  our  data  and  Eq.  (Ill- 9c).  The  resultant  value  from  this 

method  is  y  =  1.285  with  a  Curie  temperature  T  =  69.594 °K.  Most  importantly,  all  the  theo- 
- 1  C 

retical  values  of  x  for  an  exact  fit  with  these  values  fall  within  the  experimental  error  of  the 

°  _  i 

values  measured.  The  nature  of  the  agreement  is  indicated  by  the  logarithmic  plot  of  xq  vs 
Tc  —  T  shown  in  Fig.  Ill- 9.  On  the  basis  of  these  results,  we  can  conservatively  bracket  the 
value  of  the  y  by  setting  y  =  1.282  ±  0.008. 

17 

These  results  agree  with  the  scaling  law  prediction  y  =  0(6  —  1 )  to  well  within  experimental 

uncertainty.  The  scaling  laws  further  predict  that  a  =  2  —  0(6  +  1)  =  —  0.009  =*  0,  where  a  is 

defined  in  terms  of  the  magnetic  specific  heat  near  T^  by  a  |T  —  T  |”°\  For  a  =  0,  a 

logarithmic  dependence  In  |T  —  Tc  |  is  expected,  and  is  observed18  in  EuO  for  T  >  T^. 

However,  this  apparent  agreement  is  clouded  by  the  fact  that  a  logarithmic  dependence  is  not 

observed  for  T  <  T  . 

c 

An  investigation  has  been  initiated  to  determine  whether  our  data  above  Tc  can  be  fitted  to 

a  single  magnetic  equation  of  state.  Preliminary  results  indicate  that  this  can  be  done,  and  that 

1 9 

a  further  transformation  of  the  data  by  the  method  introduced  by  Schofield  yields  a  straight 
line  relationship.  Data  taken  for  T  <  Tc  is  now  being  analyzed  to  see  if  a  single  linear  relation¬ 
ship  can  describe  the  data  both  above  and  below  T  . 

c  N.  Menyuk 

K.  Dwight 


2.  Impurity  Magnon  Modes  in  KNiFgiMn 

We  have  continued  our  study  of  two- magnon  Raman  scattering  in  KNiF~  (Ref.  20)  by  using 
2+ 

samples  doped  with  Mn  impurities.  In  addition  to  the  intrinsic  two- magnon  excitation  pre- 

o 

viously  seen  in  the  pure  samples,  we  observe  a  relatively  narrow  line  shifted  from  the  5 145- A 
excitation  by  an  energy  ~550  cm"  .  This  is  well  below  the  maximum  energy  of  the  intrinsic  two- 
magnon  line,  8  50  cm  *.  Only  diagonal  elements  of  the  scattering  tensor  have  been  observed  for 
the  impurity ’excitation.  Samples  having  two  different  impurity  concentrations,  3.7  percent  and 
~2  percent  (atomic),  were  studied;  the  impurity  mode  intensity  was  stronger  in  the  higher  con¬ 
centration  sample.  This  line  disappeared  at  temperatures  approaching  T^  =  253 °K. 

From  its  temperature  dependence  and  energy,  it  seems  likely  that  this  is  a  two-magnon 

2  + 

excitation  associated  with  the  Mn  impurities,  most  probably  with  one  spin  excitation  on  the 

impurity  and  another  on  the  surrounding  intrinsic  spins.  This  situation  differs  from  the  pre- 
21 

viously  reported  case  of  KMnF^Ni  since  the  two-magnon  sq  +  d  impurity  mode  observed  there 
was  a  true  local  mode,  lying  above  the  intrinsic  two-magnon  band.  The  present  situation  in 

which  the  impurity- host  exchange  constant  is  less  than  the  host  constant,  and  in  which  an  impurity 

22 

spin  excitation  lies  within  the  host  magnon  band,  should  result  in  resonant  or  virtual  modes. 

A  complete  theoretical  treatment  of  this  problem  must  include  not  only  the  different  spin  and 

23 

exchange  of  the  impurity,  but  magnon- magnon  binding  as  well.  Since  the  approximations  which 
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21 

entered  the  theoretical  Green's  function  treatment  of  KMnF^iNi  may  not  be  valid  for  KNiF^Mn, 
we  are  presently  examining  this  theory  in  the  general  case  where  single  magnon  local  modes  may 


not  occur. 


S.  R.  Chinn 
H.  J.  Zeiger 
J.  R.  O'Connor 


3.  Raman  Scattering  in  CsMnF^ 

o 

The  Raman  spectrum  of  the  antiferromagnet  CsMnF^  has  been  observed  using  488  0- A  radia¬ 
tion  from  an  argon  ion  laser.  This  compound,  which  is  isomorphic  with  the  ferrimagnet  RbNiFv 

24  25  4  ^ 

whose  Raman  spectrum  has  also  been  observed,  '  belongs  to  the  space  group  and  has 

30  atoms  per  unit  cell.  Of  the  possible  19  Raman  active  phonon  modes2'’  we  have  so  far  tenta- 

- 1  - 1 

tively  identified  5  E^,  5  A^,  and  7  modes  ranging  in  energy  from  47  cm  to  470cm 

Several  of  these  modes  have  a  marked  temperature  dependence,  with  observed  linewidths  in- 

- 1  - 1 

creasing  from  the  instrumental  linewidth  (-3  cm  )  to  -20cm  ,  as  the  temperature  is  raised 

from  5°K  to  300 °K.  Small  energy  shifts  also  occur  in  some  of  these  modes  in  this  temperature 
range.  These  features  suggest  that  anharmonic  interactions  are  quite  large  in  this  material, 

and  we  are  presently  investigating  this  question  in  a  more  quantitative  fashion. 

2  6 

CsMnF^  orders  antiferromagnetically  below  53.5  °K.  A  Raman  mode  has  been  observed 

whose  shape,  position,  temperature  dependence  and  polarization  indicate  that  it  arises  from  the 

_  \ 

excitation  of  two  magnons.  Maximum  intensity  occurs  near  an  energy  shift  of  9  3  cm  at  low 
temperature,  and  the  line  broadens,  shifts  to  lower  energy,  and  decreases  in  intensity  as  the 
temperature  approaches  T^.  The  line  shape  suggests  that  magnon- magnon  interactions  are  also 
important  in  this  material.  The  low  temperature  energy  of  this  excitation  is  also  consistent 
with  the  approximate  energy  of  two  high- density-of- states,  zone- edge  magnons,  as  determined 


using  exchange  constants  found  by  Seavey. 
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S.  R.  Chinn 


C.  LASER  SCATTERING 


1.  Light  Scattering  from  Acoustic  Plasma  Waves  and  Single -Electron  Excitations 
in  Magnetoplasmas 


In  an  effective  mass  approximation,  the  differential  cross  section  for  the  Stokes  scattering 
of  light  from  low  temperature  electrons  in  an  n-type  spherical  band  semiconductor,  with  or  with¬ 
out  a  DC  magnetic  field  B  ,  is  given  approximately  by 


,2 

-d_2_  =.  (7 
dfidw  ,eF 


Zv^ 

4,2e2 


Im 


e(o),  q). 


(Ill-  10) 


where  and  eT  are  the  incident  and  scattered  light  polarization  vectors,  a;  =  toT— 

r  ^  i  £  2  2  I  r 

q  =  ?!  -  qF.  =  (e  An*  c  )  is  the  effective  Thomson  cross  section,  and  e(a;,  qT)  is  the  longi¬ 
tudinal  dielectric  constant.  The  light  effectively  couples  to  longitudinal  density  fluctuations  of 
the  Fermi  gas,  and  the  scattering  shows  the  peaks  of  Im  (l/e).  For  long  wavelengths  and  zero 
magnetic  field  the  cross  section  exhibits  (1)  broad  quasi- elastic  scattering  of  single  electrons 
from  just  inside  to  just  outside  the  Fermi  surface  in  the  frequency  range  0  <  a;  <  qvp,  and  (2) 
plasmon  scattering  near  cj  =  the  electron  plasma  frequency.  Both  of  these  scattering 
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processes  have  been  observed.  '  If  a  strong  magnetic  field  is  applied  to  the  electron  plasma, 
so  that  the  Landau  orbital  quantization  is  important  and  the  geometry  ~q  1 1  B0  is  considered,  the 
scattering  spectra  peaks  occur  in  the  same  frequency  bands.  One  expects  plasmon  scattering 
at  <J  =  ojp  and  quasi-elastic  scattering  in  the  range  0  <  <  qv^.  However,  the  nature  of  the 

quasi- elastic  scattering  is  drastically  altered  (see  Ref.  30  for  a  discussion  of  the  qualitative 
features  of  this  scattering).  The  electrons  on  each  occupied  Landau  level  effectively  become  a 
one-dimensional  electron  gas,  and  the  single-electron  scattering  has  multiple  peaks  at  oj  =  qvFn<j’ 
where  Vpn<j  is  the  velocity  along  Bq  of  an  electron  at  the  Fermi  surface  with  Landau  orbital 
quantum  number  n  and  spin  a.  Furthermore,  for  long  relaxation  times  the  quasi- elastic  scat¬ 
tering  will  have  additional  peaks  due  to  the  excitation  of  longitudinal  acoustic  plasma  waves  at 

frequencies  in  between  the  single  electron  peaks.  The  acoustic  plasma  wave  scattering  results 
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from  a  low  frequency  zero  of  the  longitudinal  dielectric  constant.  Several  authors  have 

independently  suggested  the  existence  of  such  waves,  but  they  have  not  been  observed.  Neither 

of  the  above  features  of  magnetoplasma  quasi- elastic  light  scattering  have  been  observed.  Al- 
34 

though  Wolff  has  given  some  theoretical  discussion  of  this  scattering  in  addition  to  that  given 
in  Ref.  30,  no  detailed  spectra  have  been  calculated  and  a  number  of  features  of  the  magneto¬ 
plasma  scattering  have  not  been  discussed.  With  these  matters  in  mind,  we  have  made  some 
preliminary  calculations  of  the  quasi- elastic  scattering  spectra  using  Eq.  (Ill- 10)  to  treat  scatter¬ 
ing  from  electrons  in  n-type  InSb. 

Because  of  the  multiplicity  of  scattering  resonances  in  the  frequency  range  [O.qVp]  the  in¬ 
clusion  of  relaxation  effects  in  the  calculation  of  c(cj,  q  )  is  important.  For  very  short  relaxation 
times  the  multiple  peaks  may  strongly  overlap,  yielding  a  broad,  smooth  spectrum  similar  to 
that  found  for  zero  magnetic  field.  For  simplicity,  we  use  a  relaxation  time  ansatz  with  a  con¬ 
stant  relaxation  time  t  in  the  equation  of  motion  for  the  density  matrix  (which  is  used  to  derive 
e).  The  use  of  the  relaxation  time  ansatz  in  quantum  transport  theory  has  some  hazards.  Very 
simple  forms  fail  (1)  to  be  gauge  invariant  and/or  (2)  to  conserve  charge  density.^  Green,  et  al.?^ 
have  given  one  reasonably  satisfactory  formulation  for  the  conductivity  problem.  However,  their 
ansatz  does  not  conserve  the  spin  density,  as  would  be  expected  for  orbital  relaxation  processes. 
We  have  modified  their  treatment  to  eliminate  this  difficulty  and  find  for  t, 
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with  lo'.cr^the  eigenstate  of  the  electron  effective  mass  Hamiltonian  [H|(*o)>=  Eo>oJa'(j>]  and 

f  the  Fermi  occupation  factor. 
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Figures  III- 10  and  III- 11  show  relative  scattering  spectra  computed  for  InSb  with  the 
following  parameters:  X1  =  10  |i,  m*/m  =  0.015,  g*  =  50,  =  16,  n0  =  1017cm“3.  The  curve 

in  Fig.  Ill- 10  was  calculated  for  a  very  long  relaxation  time,  qv^T  =  500,  and  illustrates  the 
features  described  above.  The  three  strong  peaks  each  result  from  excitation  of  an  acoustic 
plasma  wave  and  dominate  the  scattering.  The  four  small  peaks  result  from  single- electron 
excitations  on  each  of  the  four  occupied  Landau  levels.  Note  the  correlation  of  the  latter  peaks 
with  the  values  of  qvpno.  which  are  marked  by  the  arrows.  Figure  III- 11  shows  the  development 
of  the  spectra  as  the  magnetic  field  is  varied  for  a  fixed  qvpT  =10.  For  a  field  such  that  there 
are  N  occupied  Landau  levels  there  are  (N  —  1)  acoustic  waves.  Since  the  acoustic  wave  scatter¬ 
ing  is  dominant,  there  are  (N  -  1 )  scattering  peaks.  As  the  field  increases  there  are  fewer  and 
fewer  levels  occupied  and,  therefore,  fewer  scattering  peaks.  The  existence  of  an  acoustic  wave 
requires  the  occupation  of  more  than  one  level,  since,  roughly,  the  wave  involves  the  motion  of 
electrons  on  different  Landau  levels  in  a  manner  yielding  an  almost  charge- neutral  wave.  For 
fta>c/EF  =  1*5»  onlY  one  Landau  level  is  occupied,  no  acoustic  waves  exist,  and  the  scattering  is 
dramatically  reduced  to  a  broad  smear. 

Although  this  model  calculation  illustrates  some  of  the  general  features  of  magnetoplasma 
quasi- elastic  scattering,  it  is  not  expected  to  be  quantitatively  accurate  for  InSb.  A  value  of 
qvpT  =  10  is  optimistic  by  a  factor  of  five  for  10  p  light  and  typical  n-InSb.  For  qvpT  =  2,  the 
spectra  of  Fig.  Ill- 11  are  strongly  smeared  and  tend  to  become  lost  in  the  wings  of  the  plasmon 
scattering.  Also,  coupling  mechanisms  involving  the  spin  degree  of  freedom  and  spin-orbit 
coupling  are  important  in  real  semiconductors,29' 37' 30' 34  and  are  not  treated  here.  These 
mechanisms  can  significantly  alter  the  order  of  magnitude  of  the  cross  section,  the  polarization 
selection  rules,  and  the  nature  of  the  lineshapes.  A  more  realistic  calculation  extending  the 
treatments  of  Refs.  30  and  34  is  in  progress. 


F.  A.  Blum 
R.  W.  Davies 
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Fig.  Ill-l  1.  Relative  cross  section  vs  frequency  for  InSb  with  magnetic  field 
as  a  parameter  and  qVpT  =  10. 
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2.  Millimeter  and  Far  Infrared  Frequency  Mixing  in  GaAs 

Heterodyne  experiments  with  both  millimeter  and  infrared  sources  have  been  performed  in 
high  purity,  epitaxial  GaAs.  The  millimeter  waves  produce  avalanche  breakdown,  and  mixing 
via  hot  electrons  is  observed.  At  laser  frequencies  the  DC  bias  level  determines  whether  detec¬ 
tion  proceeds  via  a  similar  hot  electron  mechanism  or  through  photothermal  ionization  of  shallow 
donor  states.  The  time  constant  associated  with  the  latter  process  is  significantly  faster  than 
the  hot  carrier  relaxation. 

At  4  mm,  which  corresponds  to  ~0.3  me V,  the  response  of  GaAs  cannot  take  place  by  the 

photothermal  ionization  of  shallow  impurity  levels.  However,  since  the  conduction  electron 

- 12 

scattering  time  is  approximately  10  sec  ('—170  GHz),  the  comparatively  slow  4-mm  radiation 

can  transfer  energy  to  the  electrons  and  induce  impact  ionization.  Mixing  (and  detection)  can  then 
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proceed  via  a  hot  electron  mechanism  which  has  been  extensively  investigated  in  InSb.  Two 
4-mm  sources,  a  DX-151  klystron  and  a  CSF  carcinotron,  were  combined  in  a  4-mm  waveguide 
directional  coupler  and  the  resulting  radiation  transmitted  down  a  lightpipe  to  the  detector  im¬ 
mersed  in  liquid  helium  at  4.2  °K.  The  IF  signal  resulting  from  the  mixing  was  fed  into  a  spec¬ 
trum  analyzer  and  measured  between  0.1  and  40  MHz.  Stabilization  to  a  particular  IF  frequency 
was  provided  by  a  10- MHz  syncriminator  which  controlled  the  reflector  voltage  of  the  klystron. 
The  10  MHz  was  obtained  by  mixing  the  IF  signal  with  a  local  oscillator  in  a  balanced  diode  mixer. 
In  order  to  insure  independence  of  possible  amplitude  variations  in  the  primary  beams,  a  con¬ 
stant  IF  signal  was  maintained  on  a  reference  silicon  crystal  diode. 

A  second  method  used  to  measure  the  micro- 
Tw-i- liMiwI  wave  frequency  response  of  GaAs  was  to  sinu¬ 

soidally  amplitude  modulate  the  carcinotron  be¬ 
tween  0.1  and  2.5  MHz  by  varying  the  anode 
voltage.  This  resulted  in  a  depth  of  modulation 
of  about  10  percent.  In  order  to  obtain  response 
time  curves  for  these  high  impedance  devices, 
care  was  taken  in  all  measurements  to  load  the 
detectors  with  a  low  impedance,  wide  band 
amplifier. 

The  IF  signal  frequency  roll-off  of  InSb  and 
GaAs  at  two  bias  levels  is  shown  in  Fig.  111-12. 
Values  determined  for  InSb  are  in  good  agree¬ 
ment  with  previous  measurements  and  quite  sim¬ 
ilar  to  GaAs  at  high  bias  levels.  At  very  low 
bias  levels  a  definite  decrease  in  GaAs  time  con¬ 
stant  conforms  with  the  usual  model  that  at  high 
electric  fields  hot  carrier  effects  can  increase 
the  lifetime.  Additional  insight  into  the  origin 
of  the  carriers  active  in  microwave  mixing  is 
obtained  from  the  I-V  curves  of  Fig.  Ill- 13(a)  and 
13(b).  In  these  curves  a  dramatic  reduction  of 
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Fig.  111-12.  IF  power  versus  IF  frequency  for  InSb 
and  GaAs  at  two  DC  bias  levels. 


the  DC  bias  voltage  at  breakdown  occurs  with  in- 
creasing  microwave  power.  Since  1  mW/cm  of 
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Fig.  111-13.  I  vs  V  curves  of  GaAs.  (a)  From  left  to  right,  4-mm  radiation  of  10,  5,  and  2.5  mW 
and  no  radiation,  (b)  Same  as  (a)  with  vertical  scale  changed,  (c)  I  vs  V  curve  with  ~3  mW 
HCN  laser  radiation  (20  pA/div).  (d)  Same  as  (c)  but  with  2  pA/div.  Note  that  laser  radiation 
causes  increase  in  photoconductivity  with  little  change  in  breakdown  voltage. 


4-mm  radiation  is  equivalent  to  approximately  1  volt/cm,  the  total  electric  field  (microwave 
plus  DC)  required  for  breakdown  remains  constant.  The  lOmW  of  microwave  power  actually 
used  in  these  mixing  experiments  was  therefore  sufficient  (even  at  the  lowest  DC  bias)  to  induce 
impact  ionization  and  generate  free  carriers  in  the  detector. 

In  Figs.  Ill- 13(c)  and  13(d)  the  effects  of  HCN  laser  radiation  on  the  I- V  curves  of  GaAs  dif¬ 
fer  markedly  from  those  of  the  microwaves.  The  890-GHz  laser  radiation  cannot  cause  impact 
ionization  but  rather  a  change  in  conductivity  which  results  from  the  photo- ionization  of  shallow 
donor  levels  into  the  conduction  band.  To  further  investigate  the  response  at  these  wavelengths 
the  outputs  of  two  HCN  lasers  were  combined  in  a  TPX  beam  splitter  and  directed  onto  the  GaAs. 
The  mixed  IF  signal  was  displayed  on  a  spectrum  analyzer  as  the  lasers  were  tuned  apart  in 
frequency.  Because  of  the  limited  width  of  the  spontaneous  HCN  emission  line,  laser  oscillation 
could  only  be  maintained  up  to  a  separation  of  4  to  5  MHz.  In  Fig.  Ill- 14  the  laser  beats  are 
displayed  as  a  function  of  frequency  up  to  2.5  MHz.  The  InSb  and  GaAs  biased  above  breakdown 
show  an  IF  frequency  fall-off  below  1  MHz.  At  high  bias  the  GaAs  response  time  is  associated 
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-8-8658 


Fig.  111-14.  Output  of  spectrum  analyzer  (0  to  2.5  MHz).  Photographs  represent  a  superposition 
of  2-sec  exposures  at  each  of  5  frequency  settings.  Mixing  signals  of  two  HCN  lasers  are  obtained 
from  (a)  InSb  (bias  =  0. 2  ma),  (b)  GaAs  above  breakdown  (70pa),  (c)  GaAs  below  breakdown 
(1 .5  pa).  Response  in  the  last  case  extends  beyond  4  MHz. 

with  hot  electron  relaxation  to  the  equilibrium  distributed  in  the  conduction  band.  Below  break¬ 
down,  the  response  takes  place  via  a  photo- ionization  mechanism.  In  this  case,  after  correcting 
for  variation  in  primary  laser  power  drop-off,  the  IF  signal  shows  essentially  flat  response 
well  out  to  4.0MHz  (r  <  40 nsec).  The  recombination  to  shallow  donor  states  in  this  case  had  al¬ 
ready  been  shown  to  be  faster  than  1  psec  by  means  of  a  germanium  impact  ionization  modulator, 

40 

and  estimated  to  be  as  short  as  10  nsec  by  pulsed  ionization  measurements  in  GaAs. 

H.  R.  Fetterman 
P.  E.  Tannenwald 
C.  D.  Parker 

3.  Nonlinear  Effects  on  the  Polarization  and  Linewidth  of  Intense 
Broadband  Maser  Signals 

The  circular  polarization  properties  of  the  many  galactic  sources  of  maser  emission  by 
interstellar  OH  molecules  has  eluded  simple  explanation.  Explanations  of  three  types  have  been 
proposed:  first,  that  the  maser  is  being  pumped  by  circularly  polarized  sources;  second,  that 
systematic  Doppler  shifts  and  Zeeman  splitting  along  the  propagation  path  produces  significant 
maser  gain  for  only  one  sense  of  circular  polarization;  and  third,  that  nonlinear  competition  be¬ 
tween  opposite  circularly  polarized  signals  is  occurring,  as  in  saturated  laser  amplifiers  in  the 
laboratory.  Without  cataloging  their  shortcomings,  let  it  suffice  to  say  that  all  three  explanations 
do  not  correlate  well  with  the  observations  and  are  not  free  of  theoretical  objections.  However, 
a  new  means  of  nonlinear  competition,  parametric  down- conversion  of  one  Zeeman- split  micro- 
wave  mode  into  the  lower  frequency  microwave  mode,  and  an  electron  cyclotron  idler- wave, 
seem  to  provide  an  adequate  explanation  of  the  observations.  This  process  does  not  require,  as 
in  the  laboratory  case,  Raman-type  cross  saturation,  which  is  negligible  for  these  signals  be¬ 
cause  they  have  bandwidths  that  are  much  larger  than  the  inverse-lifetimes  of  the  maser  molecular 
states.  The  current- source  for  the  low  frequency  cyclotron  wave  is  driven  by  the  v  X  B  force  due 

to  two  phase  matching  microwave  Zeeman  modes.  Details  will  appear  soon  in  The  Physical 
41 

Review. 

Further  work  on  nonlinear  effects  deals  with  the  linewidth  of  a  broadband  maser  signal  when 

-1/2 

the  maser  amplifier  becomes  saturated.  The  expected  line  narrowing  by  a  factor  (In A)  '  , 
where  A  is  the  amplification,  gives  way  to  a  line  broadening  due  to  saturation.  The  appropriate 
spectral  power  output  for  a  saturated  two-level  maser  amplifier  has  been  calculated  from 
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approximate  solutions  to  an  integral  equation  that  is  obtained  from  the  equations  of  motion  for 
the  density  matrix  for  a  two-level  system  and  for  Gaussian  statistics  for  the  signal.  The  follow¬ 
ing  were  considered:  the  signals  are  either  quasi- monochromatic  or  broadband  relative  to  the 
inverse- lifetimes  of  the  states,  the  molecular  line  broadening  is  either  homogeneous  or  inhomo¬ 
geneous,  and  the  saturation  rate  is  either  much  less  or  much  greater  than  the  signal  bandwidth. 
The  results  for  the  case  of  inhomogeneous  broadening,  broadband  signals,  and  moderate  satura¬ 
tion  have  been  used  to  obtain  the  saturated  linewidth  and  intensity  for  the  OH  and  H^O  travelling- 
wave  masers  as  a  function  of  amplifier  length.  This  should  be  helpful  in  understanding  the 
many  observations.  M.  M.  Litvak 
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IV.  MICROELECTRONICS 


A.  ALUMINUM  BEAM  LEADS  AND  BEAM-LEAD  CROSSOVERS  ON  CERAMIC 

Beam  leads  and  beam-lead  crossovers  have  been  fabricated  on  ceramic  substrates  using 
1-3 

gold  metallization  systems.  However,  because  of  the  potential  problems  and  failure  modes 
associated  with  bonds  between  gold  beams  and  the  aluminum  metallization  on  an  integrated  cir¬ 
cuit  chip,  a  program  to  develop  a  technique  for  forming  interconnections,  crossovers  and  beam 
leads  on  ceramic  substrates  using  an  all  aluminum  metallization  system  has  been  pursued. 

The  ceramic  substrates  are  prepared  either  by  ultrasonically  drilling  the  required  holes  in 
as-fired,  99.5  percent  A12C>3  with  an  8-micro-inch  or  better  surface  finish,  or  by  having  the 
ceramic  manufacturer  prepunch  the  holes  in  green  ceramic  tape  before  firing.  The  latter  method 
is  preferred  as  it  greatly  reduces  the  amount  of  cracking  encountered  during  subsequent  process¬ 
ing  of  the  substrate. 

In  order  to  form  the  beam  leads  over  the  holes  in  the  ceramic,  the  hole  must  be  filled  or 
covered  (tented)  with  a  material  that  can  be  leached  away  without  affecting  the  aluminum.  The 
low-melting  glass  used  as  a  filler  material  for  gold  beam  lead  fabrication  cannot  be  used  because 
the  HF  used  to  leach  the  glass  will  rapidly  attack  aluminum.  A  promising  material  that  has  been 
used  to  tent  the  holes  is  Riston,  a  sheet  photoresist  product  of  DuPont.  The  photoresist  sheet 
is  applied  to  the  entire  substrate  and  exposed,  leaving  squares  slightly  larger  than,  and  covering 
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the  holes.  Using  a  filament  evaporator,  7  5,000  to  100,000  A  of  aluminum  are  deposited  on  the 
substrate.  The  interconnection  pattern  (minus  crossovers)  and  beam  leads  are  delineated  and 
etched  using  standard  photolithographical  techniques.  The  Riston  is  removed  using  either  com¬ 
mercial  photoresist  strippers  that  do  not  attack  aluminum,  or  an  oxygen-plasma  photoresist 
asher.  The  latter  method  is  much  preferred  because  of  its  speed  and  the  completeness  of  the 
Riston  removal,  a  prerequisite  for  reliable  beam  to  chip-pad  bonding. 

The  technique  used  to  apply  the  Riston  is  quite  critical.  The  ceramic  is  degreased  in  boil¬ 
ing  trichloroethylene,  rinsed  in  three  separate  beakers  of  electronic  grade  methanol,  and  blown 
dry  with  nitrogen.  The  substrate  is  then  preheated  for  15  minutes  at  80 °C  and  the  Riston  pre¬ 
heated  at  80 °C  for  2  minutes.  After  pressing  the  heated  Riston  onto  the  heated  substrate,  the 
substrate  is  placed  Riston  side  down  on  a  smooth  ceramic  plate  and  baked  for  15  minutes  at 
120°C,  and  then  baked  another  15  minutes  at  120°C  with  the  Riston  side  up.  The  substrate  is 
then  removed  from  the  oven  and  allowed  to  cool.  After  exposure  with  collimated  light  from  a 
high  pressure  mercury  arc  lamp,  the  Riston  is  spray  developed  and  the  excess  developer  blown 
off  with  nitrogen.  The  exposed  and  developed  Riston  is  post-baked  for  10  minutes  at  180°C  be¬ 
fore  loading  into  the  evaporator  for  aluminum  deposition. 

The  technique  for  forming  air-gap,  beam-leaded  crossovers  with  aluminum  is  similar  in 
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philosophy  to  that  used  for  beam-leaded  crossovers  with  gold.  However,  photoresist  KTFR  is 

used  as  the  leachable  base  for  forming  the  crossover  instead  of  copper.  The  entire  interconnec¬ 
ts 

tion  and  beam-lead  pattern  is  coated  with  40,000  to  60,000  A  of  KTFR,  and  holes  down  to  the  pads 
of  the  interconnection  pattern  where  the  crossovers  will  make  contact  are  exposed  and  developed 
in  the  photoresist.  The  substrate  is  then  baked  at  180°C  for  30  minutes  to  thoroughly  remove 

o 

all  the  solvents  from  the  KTFR.  50,000  to  7  5,000  A  of  aluminum  are  evaporated  over  the  entire 
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Fig.  IV-1.  Aluminum  beam-leaded  crossover. 


substrate,  keeping  the  substrate  temperature  at  approximately  100 °C.  The  crossover  pattern  is 
delineated  and  etched  using  standard  photolithographic  techniques,  and  all  the  photoresist  removed 
using  commercial  strippers  or  an  oxygen-plasma,  dry  photoresist  asher.  The  final  step  is  to 
sinter  the  aluminum  to  aluminum  contacts  between  the  crossovers  and  the  interconnection  pattern 
for  10  minutes  in  nitrogen  at  57 5° C.  Figure  IV-1  shows  the  aluminum  beam-leaded  crossover. 

F.  J.  Bachner 

B.  PASSIVATED  MESA  TECHNIQUE 

Most  semiconductor  devices  and  monolithic  integrated  circuits  are  fabricated  using  the  planar 
process  which  by  its  nature  affects  total  protection  of  the  sensitive  diffused  junctions  during  as 
well  as  after  various  process  steps.  Mesa  devices  have  many  desirable  advantages  over  planar 
devices  and  are  used  primarily  to  achieve  special  performance  characteristics  not  possible  with 
planar  techniques.  However,  the  mesa  process  results  in  exposed  junctions  and  no  practical 
passivation  method  has  been  developed  that  fully  protects  the  junctions.  The  process  described 
here  is  a  combination  of  the  mesa  and  planar  which  incorporates  all  the  advantages  of  the  two 
processes.  This  process  is  given  the  acronym  SIMTOP  for  the  Silicon  Nitride  (Si^N^)  Masked, 
Thermally  Oxidized,  Post  Diffused  Mesa  Process.  The  process  procedure  is  basically  as  fol¬ 
lows:  A  silicon  wafer  is  deposited  with  a  thin  layer  of  Si^N^  (Ref.  4).  The  Si^N^  is  etched5  to 
form  a  desired  geometry  using  photolithography  and  the  appropriate  etchants.  The  mesa  is  then 
formed  by  etching  the  silicon.  The  resulting  mesa  will  have  smooth  crystallographic  planes  if 
the  silicon  material  is  <100}  orientation  and  is  etched  using  an  anisotropic  etch.^”®  Then  the 
silicon  is  thermally  oxidized  to  a  thickness  of  8,000  to  10,000  A.  The  Si3N4  is  removed**  and  the 
p-n  junction  formed  by  diffusion.  The  metallized  contacts  are  formed  by  evaporating  or  plat¬ 
ing  the  appropriate  metal,  Al,  Ni-Au,  Cr-Au,  etc.  Figure  IV-2  shows  the  cross  section  of  a 
SIMTOP  device  and  Fig.  IV -3  the  top  view. 
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Fig.  IV-2.  Cross  section  of  a  SIMTOP  device. 


Fig.  IV-3.  Top  view  of  a  SIMTOP  device. 
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The  post-diffusion  feature  of  this  mesa  process  results  in  devices  which  have  all  the  advan¬ 
tages  of  the  planar  process;  high  reliability,  uniform  electrical  and  physical  characteristics  and 
fabrication  simplicity.  In  addition,  the  mesa  structure  eliminates  the  curved  junction  and  lateral 
space-change  widening  typical  of  planar  structures  and  provides  lower  capacitance  and  high  break¬ 
down  voltages.  Reliability  is  enhanced  over  the  prior  mesa  process  since  the  crystallographic 
planes  of  this  mesa  structure  tend  to  reject  impurities  and  contamination  which  would  be  incor¬ 
porated  in  the  subsequent  oxidation.  Since  the  junction  formation  is  similar  to  the  planar  p-n 
junction  formation,  no  difference  in  reliability  is  anticipated.  Thermal  properties  should  be 
enhanced  when  these  devices  are  mounted  junction  side  down  since  the  entire  mesa  top  is  metal¬ 
lized  and  contact  oxide  on  the  surface  is  not  required  as  it  is  for  deposited  glass  or  thermal  oxide 
grown  on  normal  mesa  devices.  The  thermal  oxide  in  these  SIMTOP  devices  is  confined  to  the 
slopes  of  the  mesa  only. 

An  Impatt  diode  was  fabricated  with  this  process  to  illustrate  the  feasibility  of  the  process. 
Although  the  diode  design  was  not  optimized  for  X-band  it  did  oscillate  in  X-band  in  a  CW  mode 
and  more  important,  it  did  exhibit  stable  sharp  reverse  breakdown  voltage  at  low  leakage  and 
good  forward  voltage  drop  in  either  conventional  or  flip  mounted  configuration.  The  diode  ab¬ 
sorbed  and  dissipated  5  watts  CW  power  although  it  was  a  small  (2pf  at  0  volts)  area  device.  It 
is  expected  that  with  further  work  in  the  microwave  field,  more  efficient  diodes  will  be  fabricated 

which  will  also  retain  all  the  advantages  of  the  SIMTOP  process. 

10 

New  geometrical  designs  such  as  the  ring,  star,  cross,  grid  or  comb  structures  are  pos¬ 
sible  in  mesa  form  for  devices  such  as  the  PIN,  varactor,  Impatt  and  Trapatt  elements.  The 
SIMTOP  process  also  has  considerable  advantages  in  high  voltage  devices  for  high  power  and 
high  frequency  applications.  R  A>  Cohen 

C.  COMPUTER  AIDED  DESIGN  AND  DEVELOPMENT 

Work  on  the  Mannplot  and  Digitizer  mask  making  programs  has  continued  over  the  past 
quarter.  The  programs  are  continually  being  refined. 

An  on-line  version  of  the  CIRCUS  circuit  analysis  program  has  essentially  been  completed. 
This  program  is  designed  to  allow  a  circuit  designer  to  enter  circuit  descriptions  to  the  program 
via  a  time  shared  console  and  get  back  a  limited  analysis  of  his  circuit’s  behavior  interactively. 
Complete  analyses  may  then  be  submitted  for  batch  processing  using  the  data  entered  through 
the  time  shared  program.  It  is  hoped  that  this  will  increase  use  of  the  CIRCUS  program  by  mak¬ 
ing  it  easily  accessible  to  all  laboratory  designers.  Documentation  will  be  distributed  within  the 
Laboratory  and  will  also  be  available  to  all  interested  outside  parties  as  soon  as  it  is  printed. 

N.  B.  Childs 
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